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A  Mathematical  Model  of  the  Inertial 
Properties  of  a  Carrier-Backpack  System 


INTRODUCTION 

In  the  field  of  biomechanics,  researchers  use  a  variety  of  research 
techniques  to  evaluate  various  aspects  of  physical  performance.  Mathematical 
modeling  is  one  analysis  technique  which  is  available  to  the  biomechanics 
researcher  for  relating  performance  to  the  mechanical  characteristics  of 
the  system  under  investigation. 

According  to  Olinick,^  a  researcher  may  develop  conclusions  about  a 
situation  under  observation  by  conducting  controlled  experiments  and  recording 
the  results  or  by  developing  a  model  which,  when  carefully  designed,  will 
eventually  lead  to  the  same  conclusions  as  experimentation.  The  advantage 
of  mathematical  modeling  is  that  a  model  is  free  of  the  physical  limitations 
often  associated  with  experimental  research.  Selected  variables  of  interest 
can  be  isolated  and  Carefully  controlled. 

There  are  basically  two  classifications  of  mathematical  models  which  are 
commonly  used:  deterministic  and  probabilistic.  Deterministic  models  are 
those  which  are  based  on  exact  relationships  such  that,  for  any  given  input, 
the  exact  end  result  may  be  determined.  Probabilistic  models,  as  the  name 
implies,  are  based  on  the  assumption  that  an  observation  or  system  under 
investigation  can  occupy  one  of  several  different  states  with  different 
probabilities  at  any  given  moment.  Consequently,  inferential  statistics  play 
an  Important  role  in  the  development  and  evaluation  of  these  models.  In  either 
case,  the  goal  of  the  researcher  is  to  develop  the  most  accurate,  yet  simple, 
model  possible. 

In  this  project  a  deterministic  model  was  developed  to  examine  the 
Inertial  characteristics  of  a  human-backpack  system.  In  any  system,  the  accelera¬ 
tion  of  a  body  is  related  to  the  inertial  properties  of  the  body.  In  the  case 
of  linear  motion,  the  mass  of  an  object  representr  the  resistance  to  linear 
acceleration.  When  considering  angular  motion,  however,  one  must  be  concerned 
not  only  with  the  magnitude  of  the  object’s  mass,  but  also  with  how  that  mass  is 
distributed  about  some  axis  of  rotation.  The  moment  of  inertia  of  an  object 
takes  into  consideration  both  mass  and  the  distribution  of  that  mass  and  thus 
represents  the  reoistance  of  that  object  to  angular  acceleration. 

The  Inertial  properties  of  a  human-pack  system  are  quite  important  to 
the  ability  of  an  individual  to  make  rapid  movements,  either  linear  or  angular 
in  nature.  For  example,  a  soldier  in  a  combat  situation  may  be  required  to 
make  some  evasive  maneuver.  A  greater  mass  and  greater  moment  of  inertia  of 
the  soldier-pack  system  would  Inhibit  the  ability  of  the  soldier  to  make 
rapid  linear  and  angular  accelerations,  respectively. 


^Olinick,  M. ,  An  Introduction  to  Mathematical  Models  in  the  Social  and  Life 
Sciences .  Reading,  Massachusetts:  Addison  Wesley  Publishing  Company,  1978. 
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Because  a  pack  can  be  loaded  in  various  ways,  an  individual  has  some 
control  over  the  inertial  properties  of  a  pack  and  the  influence  of  the 
pack  inertial  properties  on  the  entire  carrier-pack  system.  Consequently, 
it  was  the  purpose  of  this  project  to  examine  the  inertial  characteristics 
of  the  total  carrier-backpack  system  using  a  mathematical  model.  The  system 
was  designed  so  that  both  male  and  female  models  could  be  examined  in  conjunc¬ 
tion  with  four  different  backpacks. 

The  work  performed  in  this  project  was  an  extension  of  that  done 
previously  by  Hinrichs,  Lallemant,  and  Nelson.  In  their  work  they  examined 
the  Inertial  properties  of  throe  backpacks  of  different  design  while 
manipulating  the  loading  configurations  of  the  packs.  In  discussing  the 
results,  chey  summarized  certain  inertial  properties  which  are  desirable  in 
a  backpack,  based  on  mechanical  principles.  These  included  having  the  smallest 
possible  mass,  having  the  dorso-ventral  center  of  gravity  (CG)  component  as 
close  to  the  body  as  possible,  having  the  pack  symmetrically  loaded  from  side 
to  side,  having  the  longitudinal  CG  component  as  low  as  nossible,  and  having 
the  smallest  possible  moments  of  inertia.  The  authors  were  quick  to  point  out 
that  these  desirable  characteristics  are  not  always  feasible.  With  this 
limitation  in  mind,  their  results  suggested  that  the  most  desirable  combinations 
of  the  loading  configurations  examined  would  be  to  load  the  equipment  low  and 
place  any  added  weight  on  the  sides  and/or  the  bottom  of  the  pack.  In  addition, 
any  extra  weight  should  be  placed  as  close  to  the  pack  frame  as  possible. 

PROCEDURES 

In  this  phase  of  the  work  on  load  carrying  behavior,  a  computer  model 
was  developed  to  estimate  the  inertial  properties  of  a  human-backpack  system. 

In  order  to  understand  the  methods  used  in  this  phase,  the  reader  should  have 
some  basic  understanding  of  rigid  body  dynamics.  Because  this  information  can 
be  obtained  from  any  number  of  texts  (eg.,  Synge  and  Griffith,  1942;  Greenwood, 
1970;  or  Beer  and  Johnston,  1977),  it  will  not  be  discussed  in  this  write-up. 

The  reader  is  also  referred  to  Hinrichs  et  al.  for  discussion  of  basic  mechanical 
considerations.  Their  discussion  is  quite  appropriate  for  this  project. 

Twelve  Segment  Model 

For  the  purpose  of  this  phase  of  the  project,  a  twelve  segment  model  was 
developed  to  represent  a  soldier  wearing  a  pack  and  other  common  objects 
normally  carried  in  combat.  Information  on  the  items  included  is  presented 
in  Appendix  A.  Because  the  principal  objective  of  this  phase  wa3  to  estimate 


‘'Hinrichs,  R.N.,  S.R.  Lallemant,  and  R.C.  Nelson.  An  Investigation  of  the 
Inertial  Properties  of  Backpacks  Loaded  in  Various  Conf lguratlons .  (Tech . 
Rep.  llATICK/TR-82/G23) .  Natick,  Massachusetts:  U.S.  Army  Natick  Research 
and  Development  Laboratories,  May  1982. 
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the  influence  of  the  backpack  type  and  load  on  the  inertial  characteristics 
of  the  entire  system,  the  values  of  inertial  variables  for  all  segments, 
excep.  the  backpack  and  added  load,  were  fixed.  The  computer  model  was 
developed  so  that  one  of  four  backpacks,  ALICE  LC-2,  ALICE  LC-1,  LOCO,  and 
PACKBOARD,  could  be  used  in  conjunction  with  one,  two,  or  no  added  loads.  The 
backpack  systems  are  described  in  Appendix  A.  In  addition,  the  program  allowed 
for  the  incorporation  of  either  a  male  or  female  human  body. 

In  developing  the  model,  a  common  coordinate  system  was  used  so  that  each 
segment  could  be  appropriately  positioned  relative  to  all  other  segments.  The 
origin  of  the  coordinate  system  was  fixed  at  a  point  between  the  feet  of  the 
body  at  ground  level.  The  X-axis  represented  a  dorso-ventral  axis,  the  Y-axis 
a  transverse  axis,  and  the  Z-axis  a  longitudinal  axis.  Figure  1  gives  a  general 
representation  of  these  three  axes  relative  to  a  human  body  model. 

A  description  of  each  of  the  twelve  segments  incorporated  into  the  model 
follows. 


Human  Body  (Segment  1).  The  human  body  was  treated  as  a  rigid  body  in 
this  project.  Several  sources  of  information  were  used  to  construct  both 
male  and  female  models.  The  stature,  mass,  and  necessary  body  coordinates 
were  obtained  or  derived  using  the  anthropometric  data  from  two  separate  reports: 
The  Body  Size  of  Soldiers:  U.S.  Army  Anthropometry  -  1966^  and  Anthropometry 
of  Women  of  the  U.S.  Army  -  1977.^  For  both  the  male  and  the  female  models, 

50th  percentile  data  were  used.  These  values  were  also  adjusted  for  the  inclusion 
of  the  Army-issue  utility  shirt  and  trousers.  These  anthropometric  data  were  then 
used  in  conjunction  with  the  inertial  data  reported  by  Ilanavan^  who  developed  a 
mathematical  model  of  the  human  body  in  an  attempt  to  estimate  the  inertial 
properties  of  the  body  in  rany  different  fixed  positions.  For  this  model, 
Hanavan's  position  21,  which  is  that  shown  in  Figure  1,  was  selected.  This 
body  position  seemed  to  most  closely  approximate  the  position  of  a  soldier 
carrying  a  rifle  in  front  of  his  body.  Again  using  50th  percentile  data,  the 
location  of  the  body  center  of  gravity  and  the  body  moments  of  inertia  were 
derived  from  Hanavan's  data.  Because  the  available  inertia  values  were  for  an 
individual  with  a  body  mass  of  74,2  kg  and  a  stature  of  1.755  m,  they  had  to  be 


3 

White,  R.M.  and  E.  Churchill.  The  Body  Size  of  Soldiers:  U.S.  Army 
Anthropometry  -  1966  (Tech.  Rep.  72-81-CE).  Natick,  Massachusetts:  United 
States  Army  Natick  Laboratories,  December  1971. 

L 

Churchill,  E.,  T.  Churchill,  J.T.  McConville,  and  M.  White.  Anthropometry 
of  Women  of  the  U.S.  Army  -  1977  (Tech.  Rep.  NATICK/TR-77/024) .  Natick, 
Massachusetts:  United  States  Army  Natick  Research  and  Development  Command, 

June  1977. 

^Hanavan,  E.P.  A  Mathematical  Model  of  the  Human  Body  (Tech.  Rep.  AMRL-TR-64-102) . 
Wrlgth-Patterson  Air  Force  Base,  Ohio:  Aerospace  Medical  Research  Laboratories, 
1964. 
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adjusted  for  the  male  and  the  female  models  used  in  this  study.  These  adjustments 
were  made  by  using  correction  factors  based  on  stature  and  body  mass.  These  can 
be  shown  mathematically  as  follows: 


i  ”i2's»  i 

LA  M  2  LH 

(”h  -sa 


vv 


M  *S  ‘ 
t  A  A 
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where  Iy^,  IT^,  and  are  the  adjusted  moment  of  inertia  values  for  the 
longitudinal,  transver3e,  and  dorso-ventral  axes,  respectively;  Iyj,  ^TH> 
and  Iqh  are  the  inertia  values  from  Hanavan;  M^  and  are  the  male  or 
female  model  mass  and  stature  values;  and  M^  and  Sjj  are  the  mass  and  stature 
values  from  Hanavan.  This  adjustment  is  a  minor  one  for  the  male  model  but 
is  relatively  large  for  the  female  model.  Because  of  the  lack  of  inertia 
data  for  females,  this  adjustment  for  the  female  model  represents  a  limita¬ 
tion.  It  was  felt,  however,  that  it  would  have  a  relatively  minor  effect  on 
the  results.  The  rationale  for  this  adjustment  method  was  provided  by 
Hinrichs.6 

Helmet  (Segment  2).  Many  of  the  twelve  segments  were  modeled  as  point 
masses  rather  than  rigid  bodies.  This  greatly  simplifies  the  necessary  input 
data  since  a  point  mass  has  no  momenta  or  products  of  inertia.  This  is  an 
appropriate  assumption  as  long  as  the  transfer  term  from  a  local  coordinate 
system  to  a  coordinate  system  located  at  the  total  system  center  of  mass 
dominates  the  segment's  moments  of  inertia.  This  simplifying  assumption  will  not 
lead  to  ary  significant  error  when  applied  to  those  segments  of  relatively 
small  size  and  mass.  The  helmet  was  one  such  segment  and  consequently  was 
represented  as  a  point  mass.  The  center  of  mass  location  for  the  helmet  was 
determined  using  the  reactior-board  technique.  This  center  of  mass  location 
was  then  used  to  locate  the  helmet  relative  to  the  top  of  the  head. 

Fighting  Gear  Elements  (Segments  3-6).  Each  of  the  four  fighting  gear 
elements  was  represented  as  a  point  mass.  Rather  than  determining  center  of 
mass  positions  for  the  four  elements,  each  was  considered  as  a  regular 
rectangular  shape  with  the  center  of  mass  coincident  with  the  geometric  center. 
The  four  elements  were  then  positioned  around  the  trunk  of  the  body  using 
selected  anthropometric  measures  presented  in  the  1966  and  1977  anthropo¬ 
metric  reports.  (Ref.  3  and  A).  Figure  2  shows  the  positioning  of  the  four 
elements  in  the  transverse  and  dorso-ventral  directions.  The  location  of  the 
elements  in  these  two  directions  was  based  on  measures  of  hip  breadth  (transverse) 
and  waist  depth  (dorso-ventral).  All  four  elements  were  positioned  along  the 
longitudinal  axis  at  the  level  of  the  trochanters  of  the  femurs. 


Hinrichs,  R.N.  "Principal  Axes  and  Moments  of  Inertia  of  the  Human  Body:  An 
Investigation  of  the  Stability  of  Rotary  Motions."  Unpublished  Masters 
Thesis,  University  of  Iowa,  Iowa  City,  Iowa,  1978. 
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Figure  2.  Position  of  the  four  fighting  gear  elements  along 
the  transverse  (Y)  and  dorso-ventral ] (X)  axes. 


Boots  (Segments  7-8).  Both  boots  were  modeled  as  point  masses  as  well. 
Because  the  center  of  mars  of  the  boot  was  quite  close  to  the  approximated 
center  of  mass  of  the  foot  segment,  the  coordinates  of  the  foot  center  of 
mass  were  used  to  represent  the  location  of  the  boot. 

Rifle  (Segment  9).  Since  the  rifle  is  much  larger  and  has  a  greater 
mass  than  those  segments  treated  as  point  masses,  it  was  more  appropriately 
modeled  as  a  rigid  body.  By  using  the  reaction-board  technique,  the 
location  of  the  center  cf  mass  was  determined.  Then,  by  using  an  oscillation 
technique  which  was  previously  used  by  Hinrichs  et  al.  (Ref.  2),  the  moments 
and  products  of  inertia  were  derived.  In  order  to  define  unit  vectors  along 
the  three  local  axes  of  the  rifle,  three  points  were  established  on  the 
rifle.  Figure  3  shows  the  location  of  these  points  (A,B,C)  and  the  local 
coordinate  axes.  Point  A  was  also  used  to  position  the  rifle  relative  to 
the  body  since  it  was  assumed  that  the  right  hand  gripped  the  rifle  at 
this  point.  In  this  model  the  rifle  was  positioned  horizontally  in  front 
of  the  body.  Consequently,  knowledge  of  the  location  of  the  right  hand  on 
the  rifle  determines  the  position  of  the  left  hand  on  the  rifle.  By  knowing 
the  coordinates  of  the  right  hand  in  the  fixed  coordinate  system,  the 
location  of  the  center  of  mass  of  the  rifle  was  easily  transformed  from 
local  to  system  coordinates. 
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Figure  3.  Representation  of  the  rifle  local  coordinate  system 
and  three  points  on  the  rifle  used  to  define  unit 
vectors  along  the  three  axes. 


The  nine  segments  defined  thus  far  are  those  which  were  assumed  to  be 
fixed  in  the  computer  model.  Their  inertial  properties  are  summarized  in 
Table  1. 

Backpacks  (Segment  10).  The  origin  of  the  local  coordinate  axes  for 
each  backpack  was  located  at  the  center  of  mass  of  the  backpack  with  the 
axes  oriented  parallel  to  the  fixed  coordinate  axes.  The  local  z-axis  of 
the  backpack  was  aligned  parallel  to  the  line  connecting  the  contact  points 
of  the  pack  on  the  body  at  the  shoulders  and  hip.  This  assured  the  backpacks 
would  have  the  proper  orientation  relative  to  the  body.  The  same  reaction- 
board  and  oscillation  techniques  were  used  to  estimate  center  of  mass 
location  and  the  Inertial  properties  for  each  backpack.  For  these  tests, 
each  pack  was  loaded  with  a  sleeping  bag,  mattress,  waterproof  clothes  bag, 
poncho,  socks  and  undershirt.  These  items,  excluding  the  backpack,  totalled 
9.07  kg.  Estimates  were  then  made  of  the  location  of  each  pack  relative  to 
the  shoulder  joints.  In  the  dorso-ventral  direction,  the  edge  of  each  pack 
closest  to  the  body  was  positioned  10  centimeters  from  the  X  coordinate  of 
the  shoulder.  For  the  transverse  direction,  each  pack  was  assumed  to  be 
centered  on  the  logitudinal  axis  of  the  body.  Along  the  longitudinal  axis, 
the  tops  of  the  ALICE  LC-2,  ALICE  LC-1,  and  PACKBOARD  were  positioned  at  the 
same  level  as  the  shoulders.  Only  the  top  of  the  LOCO  was  positioned 
differently.  It  was  positioned  10  centimeters  above  the  level  of  the  shoulders 
because  of  its  greater  length.  By  representing  the  local  centers  of  mass  as 
proportions  of  the  dimensions  of  the  packs  and  then  relative  to  the  shoulders, 
the  center  of  mass  locations  were  transformed  to  the  system  coordinate  axes. 
Table  2  summarizes  the  inertial  properties  of  the  four  packs. 
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Added  Loads  (Segments  11-12).  The  computer  model  allows  for  the 
possible  addition  of  one  or  two  loads  of  variable  mass  to  the  pack  in  the 
form  of  point  masses.  The  loads  can  be  placed  in  any  position  within  the 
pack  by  specifying  the  locations  as  proportions  of  the  pack’s  three  dimensions. 
In  this  way,  a  variety  of  loading  conditions  could  be  simulated. 

Loading  Configurations 

To  examine  the  influence  of  different  loading  conditions  on  the 
inertial  prope.ties  of  the  human-backpack  s’  3tem,  seven  loading  configurations 
were  simulated  for  both  a  male  and  a  female  model  and  for  each  of  four  back¬ 
packs.  With  one  exception,  the  loading  configurations  were  simulated  under 
conditions  of  one  and  of  two  added  loads.  Each  added  load  was  6.80  kg.  Thus, 
the  pack  loads  simulated  totalled  9.07,  15.87,  and  22.67  kg.  These  values 
were  chosen  because  they  are  identical  to  pack  load  weights  used  in  other 
studies  conducted  under  this  contract. These  studies  entailed  analyses 
of  the  effects  of  load  carrying  on  the  performance  capabilities  of  men  and 
women.  For  configurations  2  through  7,  the  positioning  of  the  added  loads 
were  represented  as  a  proportion  of  the  pack  dimensions.  For  the  dorso-ventral 
direction  (X),  02  represented  the  edge  of  the  pack  closer  to  the  body.  For 
the  transverse  direction  (Y),  02  would  be  on  the  side  of  the  pack  representing 
the  models  left  side.  Finally,  for  the  longitudinal  direction  (Z),  02 
represented  the  top  of  the  pack.  The  following  is  a  brief  summary  of  the 
seven  loading  configurations.  Values  in  parentheses  represent  the  proportions 
used  in  locating  the  added  loads  in  the  X,  Y  and  Z  directions,  respectively. 

1.  No  added  load.  Both  of  the  added  loads  were  set  to  zero  so  that 
only  the  basic,  9.07-kg  load  was  included  within  the  pack. 

2.  a,b.  Performance  testing  position  (PTP) .  The  added  loads  were 

positioned  such  that  they  estimated  the  positioning  of  loads 
added  to  the  packs  during  the  performance  testing  studies 
(Refs.  7,8,9).  This  position  was  near  the  top  of  the  pack 
and  close  to  the  body.  (.10,  .50,  .20). 

3.  a,b.  High  loading.  The  added  loads  were  positioned  near  the  top 

of  the  pack  and  were  centered  in  the  pack  with  respect  to  its 
dorso-ventral  (X)  and  transverse  (Y)  axes.  (.50,  .50,  .10). 


Nelson,  R.C.  and  P.E.  Martin.  Volume  I,  Effects  of  Gender  and  Load  on 
Combative  Movement  Performance  (Tech.  Rep.  NATICK/TR-82/011) .  Natick, 
Massachusetts:  US  Army  Natick  Research  and  Development  Laboratories, 
February  1982. 


Q 

Nelson,  R.C.  and  P.E.  Martin.  Volume  II.  Effects  of  Gender,  Load,  and 
Backpack  on  Easy  Standing  and  Vertical  Jump  Performance  (Tech.  Rep. 
NATICK/TR-82/016) .  Natick,  Massachusetts:  US  Army  Natick  Research  and 
Development  Laboratories,  March  1982. 

q 

Martin,  P.E.  and  R.C.  Nelson.  Volume  III,  Effects  of  Gender,  Load,  and 
Backpack  on  the  Temporal  and  Kinematic  Characteristics  of  Walking  Gait 
(Tech.  Rep.  I«ATICK/TR-82/021) .  Natick,  Massachusetts:  US  Army  Natick 
Research  and  Development  Laboratories,  April  1982. 
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4.  a,b.  Low  loading.  The  positioning  of  the  added  loads  was  the  sane 
as  that  for  fhe  high  loading,  except  along  the  longitudinal 
(Z)  axis,  where  the  loads  were  placed  near  the  bottom  of  the 
pack.  (.50,  .50,  .90). 


5.  a,b.  Front  loading.  The  added  loads  were  centered  with  respect  to 
the  longitudinal  (Z)  and  transverse  (Y)  axes  of  the  pack,  but 
were  positioned  close  to  the  body.  (.10,  .50,  .50) 


6.  a,b.  Back  loading.  The  added  loads  were  again  centered  with  respect 
to  the  longitudinal  (Z)  and  transverse  (Y)  axes,  but  were 
positioned  near  the  edge  of  the  pack  farthest  from  the  body. 
(.90,  .50,  .50). 


7.  a,b.  Side-to-Side  loading  (S-to-S).  The  added  loads  were  posi  .ioned 
near  the  left  and  the  right  edges  of  the  pack  and  were  centered 
along  the  longitudinal  (Z)  and  dorso-ventral  (X)  axes  of  the 
body.  (.50,  .10.  .50  for  half  of  the  load  and  .50,  .90,  .50 
for  the  other  half). 


Calculated  Variables 


The  computer  model  developed  for  this  phase  of  the  project  wa3  adapted 
from  a  similar  program  written  by  Richard  N.  Hinrichs  (Ref.  6)  for  his 
investigation  of  the  stability  of  rotary  motion.  The  computer  model  generated 
values  for  the  following  variables  which  described  the  inertial  characteristics 
of  the  entire  human-pack  system:  1.  System  mass  -  the  total  mass  of  all 
components  of  the  model.  2.  System  center  of  mass  location  -  the  X,  Y,  and 
Z  values  of  the  center  of  mass  of  the  system  with  respec  to  the  fixed 
coo.'dinate  system.  3.  System  inertia  tensor  represented  at  the  system 
center  of  mass,  which  includes  the  moments  and  products  of  inertia  for  a 
ryatem  of  axes  parallel  to  the  three  fixed  axes  and  whose  origin  lies  at 
the  center  of  mass.  These  describe  how  the  total  mass  is  distributed 
about  the  three  axes.  4.  Principal  moments  of  inertia  and  their  direction 
cosines  from  the  three  fixed  axes.  These  represent  moments  of  inertia  of 
the  total  mass  about  a  new  set  of  axes  oriented  such  that  the  products  of 
inertia  are  eliminated.  The  system  inertia  tensor  and  principal  moments 
of  inertia  with  direction  cosines  provide  the  same  information  to  the 
program  user  but  in  different  formats.  For  this  reason,  only  the  values 
of  the  ine.'tia  tensor  will  be  reported  and  discussed  in  this  document. 

The  reader  is  again  referred  to  such  sources  as  Greenwood  (1970)  or  Beer  and 
Johnston  (1977)  for  details  of  inertia  tensors  and  principal  moments  of 
inertia  and  their  interpretation. 

In  addition,  the  model  was  developed  so  that  it  would  have  a  certain 
degree  of  flexibility.  The  program  user  has  the  option  of  selecting  either 
a  male  or  a  female  model  and  any  one  of  the  four  backpacks  used  in  the  testing, 
and  of  using  no,  one,  or  two  added  loads  of  variable  magnitude  positioned 
anywhere  within  the  pack.  Separate  data  decks  were  developed  for  the  male 
and  the  female  models  so  the  user  need  only  incorporate  the  appropriate  deck 
in  the  computer  analysis.  The  other  options  can  be  selected  simply  by 
adjusting  input  values  of  selected  variables  on  a  few  cards  of  the  data  deck. 
For  example t  the  subroutine  LDDAT  of  the  computer  program  reads  in  values 
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of  thM  exfra  loads  desired  and  calculates  their  location  based  on  information 
specified  by  the  user.  If  the  user  prefers  to  use  no  added  load,  the 
magnitudes  for  the  added  loads  in  the  data  deck  must  be  set  to  zero.  To 
establish  the  location  of  any  non-zero  load,  the  user  must  manipulate  the 
proportions  of  the  three  pack  dimensions  to  be  read  in  from  the  data  deck. 
These  adjustments  are  all  quite  simple  to  make  as  long  as  the  user  takes 
care  in  selecting  the  appropriate  cards  from  the  data  deck  to  be  changed. 

Although  a  certain  '.mount  of  flexibility  has  been  created  in  the 
development  of  the  computer  model,  there  are  limitations  to  the  use  of  the 
model.  This  is  particularly  true  when  considering  the  characteristics  of 
the  carrier  used  in  the  model.  The  human  body  (Segment  1  of  the  model)  was 
developed  as  being  a  50th  percentile  male  or  female  in  stature  and  body  mass. 
Because  the  emphasis  of  the  model  was  to  examine  the  effects  of  gender,  back¬ 
pack,  and  load  on  the  inertial  properties  of  the  total  system,  no  attempt 
was  made  to  incorporate  other  percentile  levels  for  the  human  body,  although 
this  could  be  done  if  desired.  In  addition,  some  limits  are  not  established 
by  the  program  itself  but  should  be  by  the  user.  For  example,  there  is  no 
specific  limit  on  the  possible  magnitude  of  added  loads  except  that  Imposed 
by  the  format  used  in  reading  this  data  in  the  program.  It  makes  little 
sense  however,  to  incorporate  loads  beyond  those  typically  used  in  the 
military. 


RESULTS  AND  DISCUSSION 


In  two  other  studies  conducted  under  this  contract,  the  effects  of 
three  main  factors  -  gender,  backpack,  and  load  -  on  selected  characteristics 
of  standing,  vertical  Jumping,  and  walking  were  examined  (Refs.  8,9).  Because 
of  the  design  of  the  computer  model  developed  for  this  study,  the  effects  of 
these  same  three  factors  on  the  interial  properties  of  the  carrier-backpack 
system  can  be  examined.  In  addition,  a  fourth  factor,  loaling  position,  is 
included.  The  influence  of  these  four  factors  will  be  discussed  In  terms  of 
three  basic  variables:  mass,  center  of  mass  location,  and  moments  of  inertia. 

Tables  3  through  10  present  results  obtained  from  executions  of  the 
computer  model.  Because  the  trends  are  quite  similar  across  the  four  factors, 
specific  results  will  not  be  discussed  for  each  condition  examined.  Rather, 
general  trends  found  in  the  data  will  be  presented. 

Gender 


Comparing  the  results  for  the  male  model  presented  in  Tables  3  to  6  with 
those  for  the  female  model  shotnn  in  Tables  7  to  10  demonstrated  that  the  male 
values  for  system  mass  and  the  three  moments  of  inertia  were  greater  than 
the  female  values  for  the  same  variables  under  all  conditions  tested.  These 
results  provide  no  new  information  since  it  is  common  knowledge  that  an 
average  male  has  a  greater  body  mass  than  an  average  female.  This  difference 
in  body  mass  is  responsible  not  only  for  the  difference  in  system  mass  but 
also  for  the  difference  in  the  values  for  moment  of  inertia. 

These  greater  values  for  the  male  model  indicated  a  greater  resistance 
to  both  linear  and  angular  accelerations  existed  for  the  male  than  for  the 
female.  This  does  not  indicate,  however,  that  It  is  easier  for  the  female  to 
accelerate  in  these  directions.  Since  the  ability  to  produce  a  linear 
acceleration  la  directly  proportional  to  the  force  acting  to  cause  the 
acceleration,  one  must  consider  the  force-producing  capabilities  for  the  male 
and  female.  An  important  and  well-known  difference  between  the  sexes  is 
the  ratio  of  strength  to  body  mass,  which  is  normally  greater  for  the  male. 
This  means  that,  even  though  the  male  has  a  greater  mass  to  accelerate,  he 
generally  has  a  greater  capacity- to  generate  force.  Consequently,  any 
adxantage  the  female  may  have  in  the  form  of  smaller  body  mass  may  well  be 
lost  due  to  a  lesser  capacity  for  generating  force.  The  same  situation 
exists  for  angular  accelerations  since  the  ability  to  produce  an  angular 
acceleration  is  directly  related  to  the  torque  generated  by  an  individual. 
Although  the  female  '  as  an  advantage  in  terms  of  a  smaller  moment  of  inertia, 
this  advantage  is  lost  due  to  a  decreased  ability  to  generate  torque. 

Finally,  a  comparison  of  the  locations  of  the  centers  of  mass  for  the 
male  and  the  female  showed  that  the  greatest  difference  existed  for  the  8 
component.  This  was  expected  since  the  coordinate  system  was  fixed  near  the 
feet  and  an  average  female  is  many  centimeters  shorter  than  an  average  male. 
Very  little  difference  existed  for  the  X  and  Y  components,  although  the  X 
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component  for  the  female  was  shifted  slightly  more  toward  the  backpack  than 
for  the  male.  This  is  due  to  the  greater  mass  of  the  pack  relative  to  total 
body  mass  for  the  female.  This  result  indicates  that  a  female  may  have  to 
make  greater  adjustments  in  posture  than  a  male  in  order  to  balance  the  pack 
and  load.  This  is  supported  by  the  results  of  the  study  of  walking  which 
showed  that  females  demonstrated  more  forward  leaning  than  males  (Ref.  9). 

Backpacks 


Again  using  comparisons  of  data  from  Tables  3  to  10,  some  general  trends 
present  for  the  packs  can  be  identified.  The  total  system  mass  is  nearly 
identical  for  all  backpacks  with  the  exception  of  the  LOCO.  Under  the 
conditions  in  which  the  LOCO  was  included,  the  total  system  mass  was  slightly 
less  than  that  for  the  other  three  packs.  This  means  that  there  is  slightly 
less  resistance  to  linear  acceleration  under  the  LOCO  conditions.  Consequently, 
in  this  sense,  it  would  be  advantageous  to  the  carrier  to  use  the  LOCO.  This 
advantage,  however,  is  quite  small. 

In  terms  of  the  moments  of  inertia,  the  results  vary  somewhat  depending 
on  the  load  and  load  position  used.  Tne  results  for  the  ALICE  LC-1  and 
ALICE  LC-2  were  nearly  identical  for  all  conditions.  This  was  expected  since 
the  packs  are  identical  and  the  frames  are  of  similar  dimensions.  In  general, 
the  LOCO  demonstrated  greater  moments  of  inertia  about  the  X  axis  than  the 
other  three  packs.  Because  of  its  greater  length,  this  was  expected.  This 
result,  however,  is  not  terribly  important  because  rotations  about  his  dorso- 
ventral  axis  are  not  regular  occurrences.  Much  more  common  are  rotations 
about  the  transverse  (e.g.  hitting  the  dirt)  and  longitudinal  (e.g.  a  change 
in  direction)  axes.  For  the  transverse  or  Y  axis,  the  LOCO  and  PACKBOARD 
appear  to  be  slightly  advantageous  over  the  ALICE  LC-1  and  ALICE  LC-2. 

Because  of  its  relatively  greater  length,  one  may  at  first  question  this 
result  for  the  LOCO  pack.  It  is  important,  however,  to  remember  that  the 
moment  of  inertia  is  dependent  upon  mass  as  well  as  the  distribution  of  the 
mass.  Even  though  the  LOCO's  greater  length  would  tend  to  add  to  moment  of 
inertia  about  the  Y  axis,  its  smaller  mass  tends  to  compensate  for  this. 

This  was  not  shown  to  be  true,  however,  for  the  X  axis.  For  the  longitudinal 
(Z)  axis,  the  LOCO  and  PACKBOARD  demonstrate  smaller  moments  of  Inertia  than 
the  ALICE  LC-1  and  ALICE  LC-2.  This  is  particularly  true  for  the  LOCO.  Not 
only  does  it  have  a  smaller  mass,  it  also  is  positioned  closer  to  the  body 
and  has  a  smaller  depth  dimension  than  the  ALICE  LC-1  and  ALICE  LC-2. 

In  terms  of  the  location  of  the  center  of  mas3,  the  LOCO  appears  to  be 
slightly  better  than  the  other  backpacks.  The  LOCO  results  in  a  smaller 
shift  in  the  direction  of  the  pack  along  the  dorso-ventral  (X)  axis  and  in  a 
lower  center  of  mass  location  along  the  longitudinal  (Z)  axis.  No  differences 
existed  between  backpacks  for  the  transverse  (Y)  axis.  The  combination  of 
these  results  suggests  that  a  slightly  more  stable  position  results  when  the 
LOCO  is  worn. I  Again,  it  is  important  to  caution  the  reader  that  all  of  these 
differences  are  rather  small  and  may  be  of  limited  practical  importance. 

Added  Loads 


The  influence  of  added  load  is  quite  simple  to  expalin.  For  a  given 
gender,  backpack,  and  load  position,  added  load  resulted  in  an  increased 
total  system  mass,  increased  moment  of  inertia  values  for  each  of  the  three 


axes,  ana  a  shift  in  the  center  of  mass  in  the  direction  of  the  added  load. 
This  is  exactly  as  was  expected  and  indicates  that  resistance  to  acceleration 
decreases  as  the  load  is  decreased.  In  other  words,  for  maximum  movement 
efficiency,  the  smallest  possible  load  should  be  carried. 

Loading  Configurations 

Unique  to  this  phase  of  the  contract  was  an  examination  of  load 
configurations.  In  the  studies  of  load  carrying  performance  of  men  and 
wome i  (Refs.  7,  8,  9),  loads  added  to  the  packs  were  added  in  only  one 
position,  (PTP).  Additional  loading  positions  were  examined  in  this  study. 
These  included  positioning  the  load  high,  low,  fr^nt  (close  to  the  body), 
back  (away  from  the  body),  and  side-to-side  (S-to-S). 

In  terms  of  the  dorso-ventral  (X)  component  of  the  center  of  mass,  the 
FTP  and  front  loadings  created  the  smallest  shift  in  the  direction  of  the 
pack.  In  both  cases  the  load  was  positioned  close  to  the  body.  As  one  might 
quickly  surmise,  the  greatest  shift  resulted  when  the  loading  was  in  the 
back  positioning.  These  results  parallel  those  for  moment  of  inertia  about 
the  longitudinal  (2)  axis.  Because  the  mass  was  distributed  closer  to  this 
axis,  the  front  and  PTP  loadings  were  best  while  the  back  loading  was  the 
poorest.  , 

Perhaps  the  most  significant  differences  between  the  loading  configura¬ 
tions  were  noted  for  the  2-component  of  the  center  of  mass  and  the  moments 
of  inertia  about,  the  dorso-ventral  (X)  and  transverse  (Y)  axes.  The  optimal 
loading  position  was  found  to  be  the  low  condition.  This  not  only  led  to  a 
more  stable  (lower)  2  component  but  also  to  smaller  moments  of  inertia  about 
the  X  and  Y  axes.  In  general,  the  poorest  loading  condition  in  terms  of  the 
2  component  and  XX  and  YY  moments  of  inertia  about  the  dorso-ventral  and 
transverse  axes  was  found  to  be  the  high  condition. 

Products  of  Inertia 

Thus  far  little  has  been  mentioned  about  the  products  of  Inertia 
resulting  from  the  different  loading  configurations.  These  results  are  not 
the  most  important  for  this  model,  but  do  provide  some  useful  information. 

In  general,  the  presence  of  a  non-zero  product  of  Inertia  indicates  that 
some  asymmetry  was  present  in  the  loading.  For  the  loading  configurations 
used  in  this  project,  the  XZ  plane  was  quite  close  to  being  a  plane  of 
symmetry  for  all  configurations.  When  this  is  the  case,  the  I^y  and  Iy£ 
products  of  inertia  are  nearly  zero.  This  situation  was  confirmed  by  the  small 


products  of  inertia  for  the  X  and  Y  axes  (I^y)  and  Y  end  Z  axes  (Iyz) , 
respectively.  The  third  product  of  inertia  (Ixz),  however,  was  considerably  J 

larger  than  the  other  two  and  represents  an  asymmetrical  loading  with  respect  i 

to  the  X  and  Z  axes.  For  example,  when  the  load  was  in  the  back  position,  t 

there  was  no  counterbalancing  load  on  the  opposite  side  of  either  the  X  or  the  S 

Z  axis.  This  resulted  in  a  significant  product  of  inertia  for  the  X  and  Z  [ 

axes.  The  larger  values  for  the  IXz  products  of  inertia  found  in  this  study  j 

confirm  this.  Even  these  values  were  relatively  small  and  do  not  indicate  ■ 

extreme  shifts  of  the  principal  axes  p*.,ay  from  the  fixed  axes  used  in  this  j 

analysis.  I 
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SUMMARY  AND  CONCLUSIONS 


In  general,  the  results  of  this  phase  of  the  project  confirmed  the 
speculations  of  Hinrichs,  Lallemant,  and  Nelson  (Ref.  2).  The  results 
demonstrated  that  the  most  desirable  loading  condition  is  one  in  which  the 
carried  mass  is  as  small  as  possible  and  is  positioned  as  low  and  close  to 
the  body  as  possible.  Loading  conditions  to  be  avoided  are  those  in  which 
the  load  is  relatively  high  and  far  from  the  body. 

In  conclusion,  computer  modeling  appears  to  be  a  very  useful  tool  for 
examining  selected  characteristics  of  lead  carrying  behavior.  The  most 
difficult  stage  in  the  modeling  process  is  the  development  of  a  simple, 
but  accurate  and  adaptable,  model.  Once  this  has  been  accomplished,  many 
different  conditions  can  be  examined  at  a  much  lower  cost  and  with  much 
less  effort  than  could  be  accomplished  through  experimentation.  This  is  not 
to  suggest  that  modeling  is  appropriate  for  all  situations.  It  can,  however, 
be  a  useful  tool  to  the  researcher. 

Recommendations  for  Further  Study 

A  logical  extension  of  this  work  would  be  examination  of  other  loading 
conditions  which  may  be  encountered  in  a  real  situation.  For  example,  it 
may  not  be  feasible  to  position  all  added  loads  in  a  symmetrical  fashion. 

If  this  is  the  case  a  probable  loading  condition  can  be  easily  simulated 
using  the  computer  model.  In  addition,  it  may  be  useful  to  develop  a  more 
refined  model  of  the  backpack  and  the  individual  components  which  commonly 
make  up  the  carried  load.  Such  a  model  would  allow  the  researcher  the 
opportunity  to  examine  many  different  loading  conditions  to  gain  greater 
insight  into  the  optimal  loading  of  a  pack.  Obviously,  this  must  be  done 
with  some  knowledge  of  what  is  needed  in  the  field.  In  other  words,  if 
modeling  suggested  that  a  particular  component  should  be  loaded  in  the  bottom 
of  the  pack  because  of  its  inertial  characteristics,  but  that  component  is 
used  quite  frequently,  it  would  be  inappropriate  to  suggest  loading  that 
component  in  any  position  other  than  an  accessible  one.  Such  a  situation, 
however,  may  well  suggest  the  need  for  a  multicomponent  pack  bag  such  that 
both  optimal  positioning  and  easy  accessibility  can  be  achieved.  Modeling 
would  be  an  ideal  tool  to  examine  such  a  question.  Finally,  inertial 
property  modeling  should  be  extended  to  include  activities  a  soldier  performs 
while  carrying  a  pack. 
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APPENDIX  A 

Clothing  and  Equipment  Used  In  This  Study 
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Clothing,  Body  Armor,  and  Sleeping  Gear 


The  items  considered  In  this  model,  either  as  components  of  the  pack 
lond  or  as  gear  worn  by  the  load  carrier,  are  standard  products  from  the 
Army's  inventory.  The  Army  nomenclature  for  each  item  and  its  military 
specification,  which  contains  a  description  of  the  item,  are  listed  below. 


Nomenclature 


Specification 


Socks,  Wool,  Cushion  Sole 

Boot,  Combat,  Leather,  Black,  Direct  Molded  Sole 
Shirt,  Utility,  Durable  Press 
Trousers,  Utility,  Durable  Press 
Undershirt,  Cotton,  White 

Helmet,  Personnel  Armor  System  Ground  Troops  (PASGT) 
Sleeping  Bag,  Intermediate  Cold,  Synthetic  Fill 
Mattress,  Pneumatic,  Insulated 
Bag,  Waterprcof,  Clothing 
Poncho,  Wet  Weather 


MIL-S-48 

MIL-B-43481E 

MIL-S-43929B 

MIL-T-43932C 

JJ-U-513D 

LP/P  DES  12-78A 

MIL-S-44016 

MIL-M-43968 

MIL-B-3108 

M1L-P-43700 


Load  Carrying  Equipment 

In  the  Army,  all  items  worn  or  carried  by  the  soldier  are  divided  into 
two  categories,  a  fighting  load  and  an  existence  load.  The  former  consists  of 
items  essential  for  the  immediate  mission,  such  as  the  clothing  and  armor  being 
worn,  a  rifle,  ammunition,  and  a  canteen.  The  existence  load  consists  of  items 
needed  to  sustain  the  soldier  in  the  field  for  a  period  of  time,  such  as 
sleeping  gear,  rations,  and  additional  clothing.  Carrying  equipment  has  been 
developed  to  accommodate  some  of  the  items  comprising  the  fighting  and  the 
existence  loads.  The  load  carrying  gear  which  was  included  in  the  present 
study  is  described  below. 

Fighting  Gear  (Figure  A-l) 

This  standard  Army  equipment  consists  of  a  belt  and  suspenders,  made  of 
nylon  webbing  and  nylon  duck,  to  which  other  items  are  attached  by  means  of  slide 
keepers.  The  equipment  hung  on  the  belt  includes: 

a.  a  cover  made  of  nylon  duck  that  holds  a  steel  cup  with  a, 9-liter 
capacity  and  a  . 9-liter  canteen  for  water. 

b.  a  plastic  case  that  holds  a  folding  intrenching  tool. 


c.  two  cases  made  of  nylon  duck  which  hold  ammunition  rounds  and  also 
have  straps  from  which  grenades  can  be  hung. 

d.  a  small  pouch  for  first  aid  dressings  or  a  compass. 

The  Army  nomenclature  and  military  specification  for  each  component  of  the 
fighting  gear  are  listed  below. 
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ALICE  Fighting  Gear 


Nomenclature 


Specification 


Belt  and  Suspenders,  All-Purpose  Lightweight 
Individual  Carrying  Equipment  (ALICE) 

Canteen,  Water,  1-Quart  Capacity 
Cup,  Water  Canteen,  Steel,  1-Quart 
Cover,  Canteen 

Intrenching  Tool,  Folding,  Lightweight 
Intrenching  Tool  Carrier 
Case,  Small  Arms,  Ammo,  30-Round 
Case,  First  Aid/Compass 

Carrying  Gear  for  Existence  Load 

Four  pack  and  frame  combinations  were  used  in  this  study.  They  included 
standard  Army,  experimental,  and  commercial  items.  Three  were  backpacks  with 
external  frames  (ALICE  LC-1,  ALICE  LC-2,  and  PACKBOARD)  and  one  was  an  internal- 
frame  system  (LOCO).  The  same  pack  was  used  on  each  of  the  external  frames. 

These  items  are  described  below. 

ALICE  Pack  (Figure  A-2).  This  standard  Army  equipment  is  a  component  of 
a  load  carrying  system  designated  as  All-Purpose  Lightweight  Individual  Carrying 
Equipment  (ALICE) .  The  ALICE  pack  is  made  of  nylon  duck  and  nylon  webbing  and 
weighs  1.3  kg.  It  has  a  large,  top-loading,  main  compartment,  an  outside  pocket 
on  each  of  two  sides  and  the  front,  and  three  smaller  pockets  above  the  center 
outside  pocket.  The  maximum  capacity  of  the  pack  is  approximately  32  kg.  The 
main  compartment  can  be  closed  by  means  of  a  drawstring  and  is  covered  by  a  storm 
flap.  The  flap  is  secured  by  two,  vertical  straps  which  encircle  the  pack.  Each 
outside  pocket  has  a  drawstring  closure  and  is  covered  by  a  flap  which  i3  secured 
by  a  single  strap.  Strips  of  webbing  sewn  on  the  outside  surface  of  the  main 
compartment  can  be  used  for  attaching  items.  A  pocket  large  enough  to  accommodate 
a  field  radio  is  sewn  inside  the  main  compartment  on  the  surface  closest  to  the 
wearer's  back.  There  are  also  "D"  rings  and  tie  strings  inside  the  main  compart¬ 
ment  which  can  be  used  to  shorten  the  pack  if  it  is  not  filled  to  capacity.  The 
pack  is  attached  to  a  frame  by  means  of  an  envelope  at  the  top  of  the  pack  which 
slides  over  the  top  of  the  frame  and  a  strap  with  a  buckle  on  the  bottom  of  each 
side  of  the  pack  which  wraps  around  the  frame. 

ALICE  LC-2  Fi.ame  (Figure  A-3) .  This  standard  Army  frame  with  its 
associated  straps  is  al3o  a  component  of  the  ALICE  system  and  is  used  with  the 
ALICE  pack.  It  carries  the  designation  "LC-2"  to  differentiate  it  from  a  frame 
(LC-1)  whic^  it  replaced  in  the  Army's  inventory.  The  ALICE  LC-2  frame  is 
structured  of  aluminum  tubing.  It  is  50.8  cm  high  and  31.1  cm  wide.  There  are 
two,  aluminum,  horizontal  members  made  from  flat  stock  which  extend  from  one  side 
of  the  frame  to  the  other  and  are  riveted  to  the  aluminum  tubing.  One,  aluminum, 
vertical  member,  also  made  from  flat  stock,  is  riveted  to  the  top  and  the  bottom  of 
the  frame.  Toward  the  top  of  the  frame,  this  vertical  piece  and  the  aluminum 
tubing  are  angled  toward  the  wearer's  back.  Two  metal  loops  are  attached  to  the 
top,  horizontal,  tubular  portion  of  the  frame.  These  are  used  to  retain  one  end 
of  the  shoulder  straps.  There  is  also  a  grommet  at  the  lower  portion  of  each 
side  of  the  frame  through  which  the  other  end  of  each  shoulder  strap  passes  and  is 
secured. 


MIL-B-43826  and 

MIL-S-43819 

MIL-C-43103 

MIL-C-43761 

MIL-C-43742 

MIL-I-43684 

MIL-I-43831 

MIL-C-4 3827 

MIL-C-43745 
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Figure  A-2.  ALICE  Pack. 
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Figure  A-2.  ALICE  Pack 


At  the  top  of  each  shoulder  strap  is  a  rectangular  piece  of  foam  spacer 
material,  22.9  cm  long,  7.0  cm  wide,  and  1.3  cm  thick,  covered  with  nylon  duck 
and  nylon  webbing.  The  remainder  of  the  strap  is  unpadded,  nylon  webbing.  A 
quick-release  device  and  a  buckle  used  for  length  adjustment  are  incorporated 
into  each  shoulder  strap.  The  lower  back  strap,  which  is  43.8  cm  long  and 
12.7  cm  high,  is  also  made  of  foam  spacer  material,  1.3  cm  thick,  covered  with 
nylon  duck.  The  back  strap  is  secured  to  the  frame  by  use  of  narrow  webbing  which 
passes  through  a  buckle.  The  waist  belt  is  comprised  of  two  pieces  of  nylon 
webbing  4.4  cm  wide.  One  end  of  each  piece  is  sewn  to  the  back  strap.  Each 
piece  includes  an  adjustment  mechanism  used  to  shorten  or  lengthen  the  belt. 

The  belt  is  secured  around  the  waist  by  a  plastic,  quick-release  device.  The 
frame  with  its  associated  straps  weighs  1.7  kg. 

ALICE  LC-1  Frame  (Figure  A-4) .  This  was  developed  for  use  with  the  ALICE 
pack  and  was  standard  Army  equipment  prior  to  the  introduction  of  the  ALICE  LC-2. 
The  LC-1  and  the  LC-2  frames  have  the  same  dimensions  and  are  of  the  same  basic 
design.  However,  the  materials  used  in  their  shoulder,  waist,  and  back  straps 
are  different.  The  top  portion  of  each  shoulder  strap,  measuring  38.7  cm  long 
aad  6.4  cm  wide,  is  made  of  a  cloth  spacer  material  covered  with  nylon  duck  . 
nylon  webbing.  The  remainder  of  the  strap  is  narrow  nylon  webbing.  A  quick- 
release  device  is  incorporated  into  the  left  shoulder  strap  and  both  straps  have 
buckles  for  length  adjustments.  The  lower  back  strap,  which  is  34.3  cm  long  ai.d 
7.6  cm  high,  is  also  made  of  a  cloth  spacer  material  covered  with  nylon  duck. 

The  back  strap  is  secured  to  the  frame  by  use  of  webbing  which  is  attached  to  a 
turnbuckle.  The  waist  belt  is  made  of  two  pieces  of  nylon  webbing  2.5  cm  wide. 

One  end  of  each  piece  is  wrapped  around  the  lower,  tubular  portion  of  the  frame. 
Each  piece  includes  a  buckle  for  adjusting  the  length  of  the  belt.  The  belt  is 
secured  around  the  waist  by  a  metal  and  plastic  quick-release  device.  The  frame 
with  its  associated  straps  weighs  1.4  kg. 

PACKBOARD  (Figure  A-5).  This  experimental  equipment,  fabricated  for  the 
study,  is  made  from  flat  aluminum  stock.  The  PACKBOARD  is  54.6  cm  high  and 
measures  34.9  cm  across  at  its  widest  point.  It  accommodates  the  ALICE  pack. 

Two  horizontal  slits  were  cut  in  the  aluminum  at  the  Cop  of  the  PACKBOARD  for 
attachment  of  the  shoulder  straps.  Two  vertical  slits  were  cut  on  each  side  toward 
the  bottom  for  attachment  of  the  lower  back  strap  and  the  straps  on  the  ALICE 
pack.  There  are  two  additional  openings  in  this  area  for  securing  the  bottom  ends 
of  the  shoulder  straps  to  the  PACKBOARD.  The  shoulder,  waist,  and  back  straps 
are  the  same  ones  used  with  the  ALICE  LC-2  frame.  A  flat,  rectangular  pad  of 
foam  spacer  material,  29.2  cm  high,  25.4  cm  wide,  and  1.3  cm  thick,  is  attached 
to  the  PACKBOARD  directly  above  the  backstrap  and  covered  with  nylon  duck. 

The  PACKBOARD  and  associated  straps  weigh  2.3  kg. 

LOCO  (Figure  A-6).  This  system  is  manufactured  by  Lowe  Alpine  Systems/ 
International  Equipment  Manufacturing.  It  is  a  top-loading,  internal-frame 
backpack.  The  frame  consists  of  two,  vertical,  aluminum  stays  which  extend  the 
length  of  the  pack,  a  distance  of  59.7  cm.  The  stays  can  be  removed  from  their 
pockets,  which  are  sewn  to  the  outside  surface  of  the  pack,  and  are  flexible  enough 
to  be  bent  by  hand.  The  stay  pockets  are  7.6  cm  apart.  The  pack  is  constructed 
of  pack  cloth  It  has  a  large  main  compartment  with  a  pocket  sewn  inside  on  the 
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surface  closest  to  the  wearer's  back.  The  :nain  compartment  can  be  closed  by 
means  of  a  drawstring  and  Is  covered  by  a  storm  flap  which  has  an  outside, 
zippered  pocket.  The  flap  is  secured  by  two  vertical  straps  and  buckles.  There 
are  three,  horizontal  straps  made  of  nylon  webbing  which  extend  aiong  each  side  of 
the  pack.  The  pack  can  be  compressed  by  use  of  buckles  attached  to  the  straps. 

A  foam  pad,  17.8  cm  high,  7.6  cm  wide,  and  .6  cm  thick,  is  attached  to  the 
center,  lower  portion  of  the  pack,  between  the  stays  of  the  frame. 

The  foam-padded  portion  of  each  shoulder  strap  is  39.4  cm  long,  6.4  cm  wide, 
and  1.3  cm  thick.  The  remainder  of  the  shoulder  strap  is  made  of  unpadded 
nylon  webbing.  The  straps  are  designed  such  that  the  padding  extends  over  the 
shoulders.  Each  strap  is  attached  to  the  pack  at  three  points.  A  strip  of 
webbing,  with  a  buckle  for  length  adjustments,  extends  from  the  middle  of  the 
padded  section  on  each  strap  to  the  top  of  the  pack.  Another  strip,  with  a 
combined  quick-release  and  length-adjustment  device,  extends  from  the  bottom 
edge  of  each  shoulder  strap's  padded  section  to  the  bottom  of  the  pack.  The 
third  attachment  point  is  at  the  center  of  the  pack,  a  location  approximating 
the  center  of  the  wearer's  back.  Here,  the  ends  of  both  shoulder  straps  are 
sewn  to  a  nylon  webbing  strap.  The  point  at  which  the  strap  attaches  to  the  pack 
can  be  adjusted  by  use  of  a  vertical  ladder  of  webbing.  A  sternum  strap  with  a 
quick-release  and  length-adjustment  buckle  extends  from  nne  shoulder  strap 
to  the  other. 

The  foam-padded  waist  belt  is  77.5  cm  long,  10.2  cm  high,  and  1.3  cm  thick. 
It  is  covered  with  pack  cloth.  Nylon  webbing  is  sewn  to  the  outside  surface  of 
the  bf’lt.  The  waist  belt  is  attached  to  the  bottom  of  the  pack  at  two  points 
(each  is  at  ‘’Upside  edge  of  a  frame  stay  pocket)  by  means  of  the  webbing  on 
the  belt,  metal  pins,  and  buckles.  The  belt  is  secured  around  the  waist  with  a 
plastic,  quick-release  device  and  webbing  straps  which  can  be  adjusted  to 
accommodate  a  range  of  waist  circumferences.  The  weight  of  the  LOCO,  including 
the  pack,  frame  stays,  and  straps,  is  1.4  kg. 

The  nomenclature  and  military  specification  for  each  pack  and  frame 
included  in  this  study  which  is  or  was  in  the  Army's  inventory  are  listed 
below. 


Nomenclature 


Specification 


Field  Pack,  Nylon,  Large,  All-Purpose  Lightweight  MIL-F-43832 

Individual  Carrying  Equipment  (ALICE) 

Straps,  Pack  Frame  and  Strap/Frame  Assembly,  MIL-S-43835 

LC-2,  All-Purpose  Lightweight  Individual 
Carrying  Equipment  (ALICE) 

Frame  Pack  with  Straps,  LC-1,  All-Purpose  MIL-F-43834 

Lightweight  Individual  Carrying  Equipment  (ALICE) 
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The  following  statement  appears  in  the  Penn  State  Computation  Center 
write-up  IMSL^®  and  contains  a  statement  regarding  the  IMSL  policies  on 
acknowledging  and  exporting  IMSL  Subroutines: 


Written  works  utilizing  computer  output  or  results  of  computational  algt-ithms 
involving  thv-  IMSL  subroutine  library  should  properly  cite  the  IMSL  Reference 
Manual.  Computer  programs  that  utilize  IMSL  subprograms  may  be  published  and 
exported  in  accordance  with  IMSL  Subscription  Policies: 

"Title  to  the  library  products  under  this  subscription 
remains  with  IMSL.  Modification  of  the  library  products 
by  the  subscriber  is  allowed,  but  such  changes  do  not 
affect  IMSL's  title  in  the  products." 

"Arrangements  with  IMSL  must  be  made  if  application 
programs  containing  embedded  object  form  IMSL  Library 
routines  are  to  be  distributed  outside  the  subscribing 
organization.  In  general,  IMSL  has  the  intent  of  being 
agreeable  in  this  case  if  the  application  of  programs 
are  to  be  made  available  at  no  charge.  If  the 
application  programs  are  to  be  made  available  for  a 
charge,  financial  arrangements  must  be  made  with  IMSL. 

When  IMSL  library  routines  are  utilized  in  research  work 
and  resulting  publications  require,  in  the  opinion  of 
the  author,  the  reproduction  of  listings  of  pertinent 
library  source  code  in  the  research  publication,  IMSL 
requests  reference,  and  grants  permission.  Where  the 
program  documentation  only  would  be  sufficient,  IMSL 
encourages  that  mode  of  usage.  When  application 
programs,  developed  as  a  by-product  of  research  work, 
contain  embedded  library  routines,  these  application 
programs  may  be  taken,  as  is,  by  the  authors  on  their 
departure  from  the  research  institution  for  internal 
usage  in  their  new  locations.  IMSL's  only  restriction 
is  that  any  embedded  library  routines  be  used  only  as 
originally  incorporated  in  the  application  program,  and 
not  removed  for  usage  in  any  other  software  development." 


1(^The  International  Mathematical  and  Statistical  Library  of  Subprograms. 

Edition  8.  University  Park,  PA:  The  Pennsylvania  State  Computation  Center, 
1980,  p.  10. 
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This  appendix  contains  a  listing  of  the  control  cards  (input  data  deck) 
for  the  male  and  the  female  models  and  the  source  listing  of  the  computer 
program.  As  one  can  quickly  see  when  examining  the  control  cards,  the  con¬ 
ditions  of  the  model  are  easily  varied  by  changing  only  a  few  cards.  Cards  1 
and  2  provide  no  information  to  vary  the  model,  but  allow  the  user  to  control 
the  number  of  different  conditions  to  be  examined  in  one  computer  run  and  to 
describe  the  various  conditions  used.  On  card  3,  the  user  must  specify 
whether  the  model  is  for  a  male  or  a  female  and  the  backpack  to  be  examined. 
This  is  the  first  of  four  cards  which  control  the  various  combinatiov.s  examined 
in  '•his  project  Cards  45,  46,  and  47  are  the  other  three  of  concern  and 
allow  the  user  to  control  the  magnitude  and  positioning  of  one  or  two  added 
loads  to  the  pack.  On  card  45,  the  user  must  specify  the  magnitude  of  the 
desired  loads  and  then  indicate  their  locations  using  cards  46  and  47.  Con¬ 
sequently,  variations  in  only  four  of  the  51  cards  of  the  control  deck  are 
responsible  for  the  many  variations  examined  in  this  project. 

When  specifying  the  magnitude  of  the  added  loads  for  a  particular  pack, 
the  user  is  liaited  only  by  the  format  statement  used  to  read  the  information 
into  the  model.  The  maximum  value  for  each  load,  therefore,  is  99.99  kg 
which  would  form  a  total  load  well  beyond  the  range  typically  used  in  the 
military.  As  far  as  the  positioning  of  the  loads  is  concerned,  the  only 
limitation  is  that  the  load  must  be  placed  somewhere  within  the  pack  since 
the  load  is  positioned  as  proportions  of  the  three  pack  dimensions.  Con¬ 
sequently,  card  46  would  contain  the  proportions  of  the  dorso-ventral,  trans¬ 
verse,  and  longitudinal  dimensions  of  the  pack  for  load  1  and  card  47  would 
contain  the  same  information  for  load  2.  One  can  see  that  the  user  has  the 
option  of  examining  almost  any  special  pack  loading.  This  illustrates  the 
major  flexibility  of  this  program. 

The  model  does  not,  however,  allow  the  user  to  examine  the  inertial 
properties  of  the  system  by  using  a  model  of  a  male  or  a  female  of  other 
than  the  50th  percentile  in  stature  and  body  mass.  Both  the  male  and  the 
female  models  were  constructed  using  50th  percentile  data  and  were  assumed 
to  be  fixed.  In  order  to  incorporate  other  body  percentiles,  one  would  need 
to  change  a  number  of  control  cards.  Values  for  stature  and  mass  (cards 
4(1),  48),  various  body  coordinate  points  (cards  17-26),  and  coordinates 
used  to  position  the  rifle  (cards  49-51)  would  all  be  affected  by  a  change 
in  the  percentile  used  for  the  male  or  the  female  model.  Given  appropriace 
anthropometric  characteristics  for  other  percentiles,  however,  the  user  could 
develop  other  male  and  female  models  and  thereby  maximize  the  flexibility  of 
the  computer  model.  Although  the  incorporation  of  different  percentiles  may 
have  some  value  in  some  future  work,  it  was  considered  to  be  beyond  the 
scope  of  this  project. 


INPUT  DATA  DECK 


Card  Variable  Values 


No. 

name 

Format 

Male 

Female 

1 

NCOND 

12 

Set  by  user 

2 

HEAD 

50A1 

Description  of 

run  supplied  by 

user 

3 

NSEX.NPACK 

211 

Both 

are  set  by 

user 

4 

SMASS 

9F6.3 

(1) 

72.655 

61.208 

(segments  1-9) 

(2) 

1.468 

1.468 

(3) 

1.198 

1.198 

(4) 

1.755 

1.755 

(5) 

1.755 

1.755 

(6) 

1.264 

1.264 

(7) 

0.852 

0.718 

(8) 

0.852 

0.718 

(9) 

3.370 

3.370 

5 

SEGMIL 

9F7.4 

(1) 

1.0710 

1.0710 

(segments  1-9) 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

r 

(4) 

0.0 

0.0 

(5) 

0.0 

0.0 

(6) 

0.0 

0.0 

(7) 

0.0 

0.0 

(8) 

0.0 

0.0 

(9) 

0.0068 

0.0068 

6 

SEGMIT 

9F7.4 

(1) 

11.7840 

11.7840 

(segments  1-9) 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

(4) 

0.0 

0.0 

(5) 

0.0 

0.0 

(6) 

0.0 

0.0 

(7) 

0.0 

0.0 

(8) 

0.0 

0.0 

'  (9) 

0.2120 

0.2120 

7 

SEGMID 

9F7.4 

(1) 

12.2450 

12.2450 

(segments  1-9) 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

(4) 

0.0 

0.0 

(5) 

0.0 

0.0 

(6) 

0.0 

0.0 

(7) 

0.0 

0.0 

(8) 

0.0 

0.0 

(9) 

0.2179 

0.2179 
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Card  Variable  yal.uea 


No. 

name 

Format 

Male 

Female 

8 

PRODIN  (Body) 

3F7.4 

'  (1) 

0,0 

0.0 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

9 

PRODIN  (Helmet) 

3F7.4 

(1) 

0.0 

0.0 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

10 

PRODIN  (FG1) 

3F7.4 

(1) 

0.0 

0.0 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

11 

PRODIN  (FG2) 

3F7.4 

(1) 

0.0 

0.0 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

12 

PRODIN (FG3) 

3F7.4 

(1) 

0.0 

0.0 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

13 

PRODIN  (FG4) 

3F7.4 

(1) 

0.0 

0.0 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

14 

PRODIN  (L.Boot) 

3F7.4 

(1) 

0.0 

0.0 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

15 

PRODIN  (R. Boot) 

3F7.4 

(1) 

0.0 

0.0 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

16 

PRODIN  (Rifle) 

3F7.4 

(1) 

0.0017 

0.0017 

(2) 

0.0041 

0.0041 

i 

(3) 

-0.0029 

-0.0029 

17 

COO  (Body) 

3F6.3 

(1) 

0.094 

0.088 

\ 

(2) 

0.0 

0.0 

1 

(3) 

0.959 

0.901 

18 

COO  (Helmet) 

3F6.3 

(1) 

0.0 

0  0 

I 

(2) 

0.0 

0.0 

1 

(3) 

1.679 

1.563 
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Card 

No. 

Variable 

name 

Format 

Male 

Values 

Female 

19 

COO(FGl) 

3F6.3 

(1) 

0.0 

0.0 

(2) 

0.205 

0.215 

(3) 

0.917 

0.836 

20 

COO  VFG2) 

3F6.3 

(1) 

0.175 

0.130 

(2) 

0.125 

0.125 

(3) 

0.917 

0.836 

21 

C00(FG3) 

3F6.3 

(1) 

0.175 

0.130 

(2) 

-0.125 

-0.125 

(3) 

0.917 

0.836 

22 

COO(FG4) 

3F6.3 

CD 

0.0 

0.0 

(2) 

-0.205 

-0.215 

(3) 

0.917 

0.836 

23 

COO (L. Boot) 

3F6.3 

CD 

0.120 

0.109 

(2) 

0.165 

0.176 

(3) 

0.0 

0.0 

24 

COO (R. Boot) 

3F6.3 

CD 

0.120 

0.109 

(2) 

-0.165 

-0.176 

(3) 

0.0 

0.0 

25 

COO (Rifle) 

3F6.3 

(1) 

0.384 

0.347 

(2) 

0.074 

0.090 

(3) 

1.067 

0.998 

26 

SHOUL 

6F6.3 

CD 

0.0 

0.0 

(2) 

0.226 

0.210 

(3) 

1.436 

1.333 

(4) 

0.0 

0.0 

C5) 

-0.226 

-0.210 

(6) 

1.436 

1.333 

27 

CF 

12F6.3 

(1) 

0.100 

0.100 

(2) 

0.0 

0.0 

(3) 

0.0 

0.0 

(4) 

0.100 

0.100 

(5) 

0.0 

0.0 

(6) 

0.0 

0.0 

(7) 

0.100 

0.100 

(8) 

0.0 

0.0 

(9) 

0.100 

0.100 

(10) 

0.100 

0.100 

(ID 

0.0 

0.0 

(12) 

0.0 

0.0 
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Card  Variable  Values 


No. 

name 

Format 

Male 

Female 

28 

PKMASS 

4F5.2 

(1) 

12.51 

12.51 

(2) 

12.09 

12.09 

(3) 

10.63 

10.63 

(4) 

12.87 

12.87 

29 

PKMI 

3F5.3 

(1) 

0.447 

0.447 

(ALICE  LC-2) 

(2) 

0.354 

•  0.354 

(3) 

0.351 

0.351 

30 

PKMI 

3F5.3 

(1) 

0.480 

0.480 

(ALICE  LC-1) 

(2) 

0.344 

0.344 

(3) 

0.345 

0.345 

31 

PKMI 

3F5.3 

(1) 

0.385 

0.385 

(LOCO) 

(2) 

0.296 

0.296 

(3) 

0.398  ■ 

0.398 

32 

PKMI 

3F5.3 

(1) 

0.540 

0.540 

(PACKBOARD) 

(2) 

0.406 

0.406 

(3) 

0.478 

0.478 

33 

PKPROD 

3F7.4 

(1) 

-0.0014 

-0.0014 

(ALICE  LC-2) 

(2) 

-0.0948 

-0.0948 

(3) 

-0.0243 

-0.0243 

34 

PKPROD 

3F7.4 

(1) 

0.0077 

0.0077 

(ALICE  LC-1) 

(2) 

-0.0812 

-0.0812 

(3) 

-0.0046 

-0.0046 

35 

PKPROD 

3F7.4 

(1) 

0.0085 

0.0085 

(LOCO) 

(2) 

0.2408 

0.2408 

(3) 

0.1368 

0.1368 

36 

PKPROD 

3F7.4 

(1) 

-0.0043 

-0.0043 

(PACKBOARD) 

(2) 

0.0558 

0.0558 

(3) 

0.0790 

0.0790 

37 

PKDIM 

3F5.3 

(1) 

0.400 

0.400 

(ALICE  LC-2) 

(2) 

0.460 

0.460 

* 

(3) 

0.510 

0.510 

38 

PKDIM 

3F5.3 

(1) 

0.390 

0.390 

(ALICE  LC-1) 

(2) 

0.460 

0.460 

(3) 

0.510 

0.510 
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Card 

No. 

Variable 

name 

Format 

Male 

Values 

Female 

39 

PKDIM 

3F5.3 

(1) 

0.280 

0.280 

(LOCO) 

(2) 

0.370 

0.370 

(3) 

0.630 

0.630 

40 

PKDIM 

3F5.3 

(1) 

0.290 

0.290 

(PACKBOARD) 

(2) 

0.460 

0.460 

(3) 

0.550 

0.550 

41 

PKPROP 

3F4.2 

(1) 

0.49 

0.49 

(ALICE  LC-2) 

(2) 

0.50 

0.50 

(3) 

0.51 

0.51 

42 

PKPROP 

3F4.2 

(1) 

0.50 

0.50 

(ALICE  LC-1) 

(2) 

0.50 

0.50 

(3) 

0.52 

0.52 

43 

PKPROP 

3F4.2 

(1) 

0.38 

0.38 

(LOCO) 

(2) 

0.50 

0.50 

(3) 

0.67 

0.67 

44 

PKPROP 

3F4.2 

(1) 

0.47 

0.47 

(PACKBOARD) 

(2) 

0.50 

0.50 

r 

(3) 

0.49 

0.49 

45 

SMASS 

2F5.2 

Set  by  user 

(added  Loads) 

• 

46 

LDPROP 

3F4.2 

Set  by  user 

(added  Load  1) 

47 

LDPROP 

3F4.2 

Set  by  user 

(added  Load  2) 

48 

STAT,  XMASS 

2F6.3 

(1) 

1.744 

1.628 

(2) 

72.655 

61.208 

49 

RIFPT1 

3F6.3 

(1) 

0.384 

0.347 

*  (2) 

-0.226 

-0.210 

(3) 

1.067 

0.998 

50 

RIFPT2 

3F6.3 

(1) 

0.384 

0.347 

(2) 

0  506 

0.552 

(3) 

1.067 

0.998 

51 

RIFPT3 

3F6.3 

(1) 

0.384 

0.347 

(2) 

-0.044 

-0.028 

(3) 

0.949 

0.880 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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THE  BIOMECHANICS  OF  LOAD  CARRYING  BEHAVIOR: 

SYSTEM  INERTIAL  CHARACTE  RISTICS  — 3D  PROGRAM. 

COMPUTER  LANGUAGE:  FORTRAN  WATPOR 

THIS  PROGRAM  HAS  ORIGINALLY  WRITTEN  BY  RICHARD  N.  HINRICHS 
FOR  dlS  MASTERS  THESIS  ENTITLED  "PRINCIPAL  AXES  AND 
MOMENTS  OF  INERTIA  OF  THE  HUMAN  BODY:  AN  INVESTIGATION 
OP  THE  STABILITY  OF  ROTARY  MOTIONS"  (UNIVERSITY  OF  IOWA, 
1978).  IT  HAS  ADAPTED  POR  USE  IN  THIS  PROJECT  BY  PHILIP 
E.  MARTIN  AND  RICHARD  N.  HINRICHS  OF  THE  PENNSYLVANIA 
STATE  UNIVERSITY  BIOMECHANICS  LABORATORY. 


THIS  PROGRAM  READS  XYZ  COORDINATES  OF  SEVERAL  BODY  AND 
EQUIPMENT  LANDMARKS  AND  COMPUTES  THE  FOLLOWING: 

A.  XYZ  COORDINATES  OF  THE  WHOLE  SYTEM  CG 

B.  THE  3  CENTRAL  MOMENTS  AND  3  CENTRAL  PRODUCTS  OF  INERTIA 
IN  XYZ  OF  THE  WHOLE  SYSTEM 

C.  THE  3  CENTRAL  PRINCIPAL  MOMENTS  OF  INERTIA  OF  THE 
WHOLE  SYSTEM 

D.  THE  ORIENTATIONS  OF  THE  3  CENTRAL  PRINCIPAL  AXES 


THE  MODEL  IS  COMPOSED  OF  12  SEGMENTS,  NUMBERED  AS  FOLLOWS 
(THE  INFORMATION  IN  PARENTHESES  INDICATES  HOW  THE  SEGMENT 
WAS  MODELED,  IE.  AS  A  RIGID  BODY  OR  A  POINT  HASS)  : 

1.  WHOLE  BODY  PLUS  CLOTHING  (RIGID) 

2.  HELMET  (POINT) 

3.  FGHT  GEAR  1  — SHOVEL  (POINT) 

4.  FGHT  GEAR  2— AMMO  PACK  (POINT) 

5.  FGHT  GEAR  3--AMMO  PACK  (POINT) 

6.  FGHT  GEAR  4--CANTEEN  (POINT) 

7.  LEFT  BOOT  (POINT) 

8.  RIGHT  BOOT  (POINT) 

9.  RIFLE  (RIGID) 

10.  PACK, FRAME,  AND  BASIC  20*  LOAD  (RIGID) 

11.  ADDED  LOAD  1  (POINT) 

12.  ADDED  LOAD  2  (POINT) 


BOTH  MALE  AND  FEMALE  MODELS  WERE  DEVELOPED  FOR  THIS  PROJECT. 
THESE  WERE  CONSTRUCTED  USING  THE  ANTHROPOMETRIC  DATA  PROVIDED 
BY  THE  1966  AND  1977  ARMY  ANTHROPOMETRIC  REPORTS.  IN  BOTH 
CASES,  THE  50TH  PERCENTILE  DATA  WERE  USED  IN  CONJUNCTION 
WITH  THE  RESULTS  OP  HANAVAN  AND  CLAUSER  ET  AL. 


THE  FOLLOWING 
CF: 


COO: 

HEAD: 

LDPROP: 

NCOND : 

NPACK : 

NSEX: 

PISCG: 


DEFINES  THE  MAJOR  VARIABLES  USED  IN  THIS  PROGRAM: 
CORRECTION  FACTORS  USED  IN  LOCATING  PACKS  RELATIVE 
TO  THE  SHOULDER. 

XYZ  COORDINATES  OF  THE  SEGMENTAL  CG'S. 

HEADER  INFORMATION  DESCRIBING  RUN. 

DESIRED  LOCATION  OF  THE  ADDED  LOADS  REPRESENTED 
AS  A  PERCENTAGE  OF  THE  PACK  DIMENSIONS. 

NUMBER  OF  DIFFERENT  CONDITIONS  TO  BE  EXAMINED 
IN  ONE  RUN. 

PACK  TO  BE  USED-- TH EY  A3  g  NUMBEBED  AS  FOLLOWS: 

1. ALICE  LC-2  2. ALICE  LC-1  3. LOCO  4 . PACKBOARD 

SEX  OF  THE  MODEL  HUMAN:  1 . MALE  2. FEMALE 
INERTIA  TENSOR  OF  EACH  SEGMENT  RELATIVE  TO  LOCAL 


58 


non 


C  AXES. 

C  PKDIH :  DIMENSIONS  OP  THE  PACKS. 

C  PKMASS:  HASS  OP  THE  PACKS. 

C  PKHI:  HOHENTS  OP  INEBTIA  OP  THE  PACKS. 

C  PKPBOD:  PRODUCTS  OP  INEBTIA  OP  THE  PACKS. 

C  PKPBOP:  CG  LOCATION  OP  EACH  PACK  REPRESENTED  AS  A  PER- 

C  CENTAGE  OP  THE  PACK  DIHENSIONS. 

C  PBODIN:  SEGMENT AL  PRODUCTS  OP  INEBTIA:  1.XT  2.XZ  3.YZ 

C  RIFPT*:  THREE  POINTS  ON  THE  RIFLE  USED  TO  DEFINE  ITS 

C  LOCATION  IN  SPACE. 

C  SEGNID:  SEGMENTAL  MOMENT  OF  INEBTIA  ABOUT  ITS  DORSO- 

C  VENTRAL  AXIS. 

C  SEGMIL:  SEGMENTAL  MOMENT  OF  INEBTIA  ABOUT  ITS  LONG- 

C  ITUDINAL  AXIS. 

C  SEGMIT:  SEGMENTAL  MOMENT  OP  INERTIA  ABOUT  ITS  TRAMS- 

C  VERSE  AXIS. 

C  SHOUL:  XIZ  COORDINATES  OP  THE  LEFT  AND  RIGHT  SHOULDERS. 

C  SHASS:  SEGMENTAL  MASSES. 

C  STAT:  BODY  HEIGHT. 

C  STSMS:  TOTAL  SYSTEM  HASS  (MOT  BODY  HASS). 

C  TISBCG:  INERTIA  TENSOR  OF  EACH  SEGMENT  RELATIVE  TO  THE 

C  SYSTEM  CG. 

C  NBCPHI:  SYSTEM  PRINCIPAL  MOMENTS  OF  INERTIA. 

C  NSICGB:  WHOLE  SYSTEM  INERTIA  TENSOR  THROUGH  SYSTEM  CG. 

C  XMASS:  HASS  OP  THE  BODY  PLUS  CLOTHING  (NOT  SYSTEM  MS). 

C  XXCG:  XYZ  COORDINATES  OP  THE  SYSTEM  CG. 

C 

C  THE  FOLLOWING  DESCRIBES  THE  UNITS  USED  ON  THE  VARIOUS 
C  TYPES  OF  VARIABLES  USED  IN  THIS  PROJECT.  THESE  UNITS 
C  APPLY  TO  BOTH  INPUT  AND  OUTPUT  DATA. 

C  DISTANCES:  METERS 

C  HASSES:  KILOGRAMS 

C  HOHENTS  OF  INERTIA:  KILOGHAM*HBTERS**2 

C  PRODUCTS  OF  INERTIA:  KILOGRAM*METERS**2 

C  PRINCIPAL  AXES:  DENOTED  BY  DIRECTION  COSINES 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
COMMON  /INEH/  SEGMIL ( 12) , SEGMIT (12) , SEGHID(12) , SHASS (12) # 

6  PRODIN  ( 12,3) 

COMMON  /COOR/  COO  (2, 1 2 ,3)  ,  SHOUL  (2,  3)  , RIFPT1  (3)  ,  HIFPT2 (3)  , 

S  HIFPT3(3) 

COMMON  /COR R P/  CP  (4,3) 

COMMON  /PK/  PKHI  (4,3)  , PKMASS  (4)  ,PKDIH  (4,3)  ,PKPROP(4,3)  , 

S  PKPROD  (4, 3) 

DIMENSION  XXCG (3) ,USL (3) ,OST (3) ,USD(3) , PISCG ( 12,3,3) ,XLAHDA (3,3) , 

S  A  (3)  ,  TISBCG  (12,3,3)  ,WfiICGR  (3,3)  ,WBCPMI(3)  ,BBC?AX(3,3)  , 

®  XL  AMDT  (3,3)  ,  PIS  CGI  (3,  3)  ,RISCG1  (3,3)  ,ZZ(3,3)  ,OP  (3) 

DOUBLE  PRECISION  XBICGR  (6)  ,  EIGNVL  (3)  ,  EIGMYR  (3  ,3)  ,WK  (  10) 

LOGICAL*1  HEAD (50) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
READ  NUMBEB  OP  CONDITIONS  TO  BE  EXAMINED  IN  THIS  RUN 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

READ (5, 10)  NCOHD 
10  FORM A 2 (12) 

C 

DO  999  NC*1 , NCOND 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 
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C  READ  U2ADE3  INFORMATION  AND  THE  SEX  AMD  PACK  TO  BE  USED 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

HEAD  (5,20)  H2AD,NSEX, HPACK 
20  FORM  AT  (50H  ,/,  211) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  MASS,  HONE 2?  MID  product  of  inertia  values,  and  CG  LOCATIONS  ABE 
C  READ  IN  P03  THE  NINE  SEGMENTS  EXCLUDING  PACK  AMD  ADDED  LOADS 

c 

C  J  INDICATES  SEGMENT  NUMBER 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

BEAD  (5,  30)  (SMASS(J)  ,J=1,9) 

30  FORMAT (9?6. 3) 

READ  (5,40)  (SEGtllL  (J)  , J=1,9)  ,  (SEGMIT  (J)  ,J  =  1,9)  ,  (SEGSID{J) 

40  FORMAT  (9F7.  4,2  (/,9F7.4)  ) 

DO  50  J=  1 , 9 

READ  (5,60)  (PRODIN  (J,K)  ,K=1,3.t 
60  FORMAT  (3?7.  4) 

50  CONTINUE 

DO  70  1=1,9 

READ  (5,  30)  (COO  (N S  EX,  J  ,  K)  , K=  1 , 3 ) 

80  FORMAT  (3F6. 3) 

70  CONTINUE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  READ  IN  COORDINATES  FOR  THE  SHOULDERS 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
READ  (5,90)  (  (SEOUL  (J,K)  ,K=  1 , 3 )  ,  J=  1 , 2) 

90  FORMAT  (6P6. 3) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  READ  IN  CORRECTION  FACTORS  USED  IN  LOCATING  PACKS  AND  ADDED 
C  LOADS  RELATIVE  TO  THE  SHOULDERS 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
READ  (5,100)  ((CF  (I,J)  ,J»1,3)  ,1=1,4) 

100  FORM  AT ( 1 2F6. 3) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  CALL  SUBROUTINES  PKDAT  AND  LDDAT  TO  READ  IN  DATA  FOR  THE  PACKS 
C  AMD  ADDED  LOADS  AND  TO  CALCULATE  LOCATIONS  FOR  THEN 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL  PKDAT ( NPACK, NS EX) 

CALL  LDDAT  (NPACK, NS  EX ) 

CCCCCO:CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  THE  STANDING  HEIGHT  (STAT)  IN  METERS  AND  THE  MASS  OF  THE  BODT 
C  ( XMASS)  IN  EG  ARE  READ  IN 

C 

CCCCCCCCCCCCCC.CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
READ  (5, 1  10)  STAT, X  MASS 
110  FORM  AT ( 2  ?6. 3) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  CORRECTION  0?  THE  BODY  MOMENTS'  OF  INERTIA: 

C  THE  UOU ENT  OF  INERTIA  VALUES  USED  ARE  FOR  AN  INDIVIDUAL  OF  73.4 
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C  KG  HASS  AMD  1.755  METERS  HEIGHT.  THESE  VALDES  SHOULD  BE  COORECTBD 
C  FOB  THE  MODEL  BEING  ANALYZED.  COHBPA  IS  THE  COORECTIOM  FACTOR  FOB 
C  THE  TRANSVERSE  AND  DORSOVE1TB1L  MOMENTS.  COBBFB  IS  THE  CORRECTION 
C  FACTOR  FOR  THE  LONGITUDINAL  MOMENT. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CORRFA»XHASS*(STAT**2) /74.2/( 1. 755**2) 

CORRFB*  (XHASS**2)  *1.755/(74.2**2)  /STAT 
SSGHIL(1)=SEGMIL(1) *CORRFB 
SEGHIT  (1) *SEGNIT (1 ) *CORRFA 
SEGHID  ( 1)*SEGNID{1)  *CORBFA 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  THE  LOCATION  OF  THE  WHOLE  SISTEM  CENTER  OF  GRAVITY  (IXCG)  IS 
C  DETERMINED 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  120  1*1.3 
XXCG (I) *0. 

120  CONTINUE 
SYSMS=0. 

DO  130  J*1, 12 

SY SB S* SYS MS* SM ASS  (J) 

130  CONTINUE 

DO  140  K* 1, 3 
X*0. 

DO  150  J*1, 12 

X*X+  SMASS  (J)  *COO  (  NSEX  ,  J  .  K) 

150  CONTINUE 

XXCG (K) =X/SYSHS 
140  CONTINUE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  READ  IN  THREE  RIFLE  COORDINATES  WHICH  DEFINE  POSITION  OF  THE 
C  RIFLE  RELATIVE  TO  THE  BODY  AND  WHICH  WILL  BE  USED  TO  DEFINE 
C  UNIT  V2CT0BS  ALONG  THE  THREE  AXES  OF  THE  RIFLE. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
READ  (5,160)  (RIFPT1  (I)  ,1*1,3)  ,  (RIFPT2  (I)  , I*  1 , 3)  ,  (RIFPT3(I)  ,I»1,3) 
160  FORM  AT ( 3F6. 3,2 (/, 3F6. 3) ) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  CALL  SUBROUTINE  AXLOC  TO  DEFINE  UNIT  VECTORS  (DSL,UST,DSD) .  THESE 
C  AXES  HERE  ESTABLISHED  DURING  OSCILLATION  TESTING. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL  AXLOC (RIFPT1,RIPPT2,RIFPT3,USL,UST,USD,0P) 

DO  170  1*1,3 
OST  (I)  =— UST  ( I) 

USD (I) *-USD (I) 

170  COMTINUB 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  THE  INERTIA  TENSOR  (PISCG)  IS  ESTABLISHED  FOR  ALL  SEGHENTS 

C  PBLATIVE  TO  THE  LOCAL  AXES  OF  EACH  SEGHENT 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  180  J*1 ,12 
PISCG  (J,  1,  1)  *SEGN ID  (J) 

PISCG  (J  ,  1  ,2)  *PRODIN  (J,  1) 
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PISCG  (J,  1,3)  =  PROD  IN  (J,  2) 

PISCG  (J,  2,  1)  =PRODIN (J,  1) 

PISCG  (J,2,2)  =SEGMIT  (J) 

PISCG  (J  ,  2, 3)  =PHODIN(J,3) 

PISCG  (o  ,3,1)  =  PRODIN(J,2) 

PISCG  (J  ,3,  2)  =PRODIN  (J,  3) 

PISCG  (J  ,3,3)  =SEGHI  L  (J) 

180  CONTIVOE 

CCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c  T«. RIFLE  INERTIA  TENSOR  BUST  BE  RELATED  TO  THE  XTZ  COORDINATE 
C  SYSTEH.  XLAMDA,  THE  TRANSFORMATION  MATRIX  BETWEEN  THE  FIXED 
C  XIZ  REFERENCE  SYSTEM  AND  THE  RIFLE  REFERENCE  SYSTEH,  AND 
C  XLAMDT,  ITS  TRANSPOSE,  ARE  DEFINED 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  190  K=1 , 3 
XLAHDA  («,  1)  =  USD  IK) 

XLAMDA  (K  ,2)  =UST  (K) 

XLAMDA  (  K,  3)  =  USL  (K) 

190  CONTINUE 

CALL  TRANSP (XLAHDA, XLAMDT) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c  THE  RIFLE  INERTIA  TENSOR  (PISCG)  IS  ROTATED  INTO  THE  FIXED  XYZ 
C  REFERENCE  FRAME 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  200  K=  1 , 3 
DO  210  L= 1 , 3 

210  PI3CG1  (K,L) =PISCG (9,K,L) 

200  CONTINUE 

CALL  HATRIH (ILAMD A , PI SCG 1 , ZZ , 3, 3, 3) 

CALL  HATBXM  (ZZ,XLAMDT,  R.ISCG1,3,3,3) 

DO  220  K*1, 3 
DO  230  L=1,3 

230  PISCG  (9  ,  K,  L)  =RISCG  1  (K  ,  L) 

220  CONTINUE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  THE  SEGMENT  INERTIA  TENSOR  (PISCG)  IS  TRANSFERED  ?IA  THE  PARALLEL 
C  AXIS  THEOREM  TO  THE  SYSTEM  CG.  TISBCG  IS  THE  SEGMENT  INERTIA 
C  TENSOR  RELATIVE  TO  THE  SYSTEM  CG  IN  XYZ 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  240  J=1 , 1 2 
DO  250  K  =  1,  3 

A  (K)  =  XXCG  (K)  -COO  ( NSEX,  J ,  K) 

250  CONTINUE 

TISBCG  (J,1,1)=PISCG(J,  1,1)  ♦SHASS  (J)  *  (A  (2)  **2+A(3)**2) 

TISBCG  (J,1, 2) =PISCG(J,  1,2)  -  SHASS(J)  *  A  (1)  *A(2) 

TISBCG  (J,1,3)=PISCG(J,  1 , 3)  -  SM  ASS  (J)  *  A  (1)  *A  (3) 

TISBCG (J, 2,1) =TISBCG(J,  1,2) 

TISBCG  (J,  2,  2)  =PTSCG  (J,  2, 2)  ♦SHASS  (J)  *  (A  (1)  **2+A(3)  ♦’♦2) 

TISBCG  (J, 2, 3)  =PI3CG  (J,  2,3)  -SH ASS  (J)  *  A  (2)  *  A ( 3) 
TISBCG(J,3,1)=TISBCG(J,  1,3) 

TISBCG  ( J,3,  2)  =TISBCG  (>3,2,3) 

TISBCG(J,3,3)=PISCS(J,3,3)  ♦  SHASS(J)*  (A(1)  ♦*2^A(2)**2) 

240  CONTINUE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
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c 

C  THE  CONTRIBUTION  OF  KACH  SEGMENT  IS  ADDBD  OP.  WBIC6H  IS  THE 
C  SYSTEM  INERTIA  TENSOR  THROOGH  ITS  CG  IN  XYZ 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

DO  260  K=1,3 
DO  270  L»1,3 
270  HBICGR  (K,L)  -0. 

260  CONTINUE 

DO  280  J*1,12 
DO  290  K*1,3 
DO  300  L=1,3 

300  NBICGR  (K,L) -WBICGB (K,L) +TISBCG ( J, K,L) 

290  CONTI NOE 
280  CONTINUE 

XBICGH(1)*WBICGR  (1,1) 

XBICGR  (2)  -NBICGR  (2,  1) 

XBICGR  (3)*WBICGR (2,2) 

XBICGR  (4)  “HBICGR  (3,  1) 

XBICGR (5) “WBICGR (3,2) 

XBICGR  (6)  *WBICGR  (3,3) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


c 

C  THE  SYSTEN  INERTIA  TENSOR  WBICGR  IS  DI1GOMALIZED  TO  DETERMINE 
C  THE  SYSTEM  PRINCIPAL  MOMENTS  OF  INERTIA  AND  THE  CORRESPONDING 

C  PRINCIPAL  AXES  DIRECTIONS.  THE  IHSL  SUBROUTINE  EIG1S  DETERMINES 

C  THE  EIGENVALUES  AND  THE  EIGENVECTORS  OP  A  REAL  SYMMETRIC  MATRIX. 

C  THE  EIGENVALUES  ARE  THE  SYSTEM  CENTRAL  PRINCIPAL  MOMENTS  OF 

C  INERTIA  (HBCPHI).  THE  EIGENVECTORS  ARB  THE  UNIT  VECTORS  IN  THE 
C  PRINCIPAL  AXES  DIRECTIONS  (NBCPAX) . 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL  SIGNS (XBICGR , 3, 2, EIGMVL, EIGN7R, 3 , WK, IER) 

DO  330  K-1,3 
DO  340  L*1,3 
WBCPAX  (K,L)  =  EIGNVE  (K,L) 

340  CONTINUE 

MBCPHT  (K)-EIGNVL(n) 

330  CONTINUE 


PERFI=WK(1) 

IERROR-IBR 

ccccccccccccccccccccccrccccccccccccccccccccccccccccccccccrcccccccccccccc 

I 

PRINT  THE  RESULTS  OF  THE  ANALYSIS 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

WRITE  (6,601) 

601  PORMAT('1', 131  (•_')) 

WRITE  (6,600)  Hi!  AD, NSEX,XMASS,STAT,SYSHS,NPACK,SHASS  (11) , 

6  SMASS  (12) 

600  FORHAT(*-*,T30, 'SYSTEM  INERTIAL  CHARACTERISTICS  FOR  A  HUMAN*, 

G  '-PACK  MODEL',///, IX, 50A1 ,//,5X, *  SEX: » , 

S  2X, 11,101,' (NOTE:  1=HALE  2-FEMALE) ' ,/, 

G  5X  ,  *  BODY  MASS  (KG)  :»  ,2X,F6. 2, /,5X,  'BODY  HEIGHT  (M)S',2X, 

G  F5. 3, /,5X, 'SYSTEM  HAS.S  (KG) : • , 2X, F6. 2./, 5X, • PACK: • , 2X,I 1, 

G  10X,  *  (NOTE:  1  =  ALICE  LC-2,  2-ALICE  LC-1,  3-LOCO,  4-PACKBOARD) ' ,/, 

G  5X, 'ADDED  LOAD  #1  MASS  (KG)  :  •  , 2X, F5. 2,/,5X, 

G  '  ADDED  LOAD  t2  HASS  (KG)  :  '  ,  2X,F5.  2) 

WRITE  (6,610) 

610  PORMAT('-' ,16X, 'SEGMENT*, 9X,»CH  COORDINATES', 10X, 
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t  ' MOMENTS  OF  INERTIA ', 121 ,' PRODUCTS  OF  INERTIA',/, 

6  5X, • SEGMENT'  ,6 X , • M ASS' , 9 X , • X • , 7 X, * Y » , 7 X, ' 2 ' ,7 X , • XX  * , 

P  8X,'YY',8X,'ZZ',8X,'XY',8X,'XZ',8X,'YZ') 

WRITE  (6,620) 

620  FORM AT (*♦',41,105 ('  •),/) 

DO  615  J=*1 , 1 2 

WRITE  (6, 630)  J,SMASS  (J)  ,  (COO  ( NS  FX  ,  J,  K)  ,K=1,3)  ,SEGMID(J)  , 

6  3EGMIT  (J)  ,  SEGMIL(J)  ,  (PRODIN(J,K)  ,K=1,3) 

630  FORHAT(8X,I2,6X,F6.  3,6X,2(P6.  3,2X)  , F6. 3 , 6 (3X , P7 . 4) ) 

615  CONTINUE 
C 

WRITE  (6, 640)  (XXCG  (I)  ,1=1,3) 

640  FORMAT  ( '  0  '  ,4  X,  '  TOTAL  SYSTEM  CM  LOCATION—  X:',P6.3,2X, 

6  '  Y:',F6.3,2X,'Z:',F6.3) 

WRITE  (6,650) 

650  FORHAT('-', 'INERTIA  T2NSOR:  MOMENTS  AND  PRODUCTS  OF  INERTIA', 

£  //,4X, 'IXI',7X,'ITT',7X,'IZZ',7X, 'IXT' ,7X, 'IXZ' ,7X, 'IIZ',/) 

WRITE  (6,660)  (K  BICGR  (K,  K)  ,F  =  1 ,  3)  ,WBICGB  (1 ,2)  ,WBICGB  (1,3)  , 

6  W  BICGR  (2 , 3) 

660  FORHAT(1X,6 (F7.3,3X))  ! 

WRITE  (6,670) 

670  FORHATC  1»,  4  (•-'),  'PRINCIPAL  MOMENTS ',  3  (•-'),  4X  ,  34  (•-•)  , 

6  'PRINCIPAL  AXES',  33('-')) 

WRITE  (6,680) 

680  FORM  AT ( ' O' ,  3X, ' 1 1 1 • , 5X , • 122 • , 5X , • 13 3* , 8X , • 1 X • , 7X , • 1 Y • , 

6  7X,  •  1Z',8X,  *2X*  ,7X,  «2Y',7X,'  27,'  ,8X  ,  •  31 '  ,  7  X,  •  3Y*  ,7X,'3Z', 

6  7  X , • PEEFI ' ,2X, ' ERR ' ) 

WRIT  E  (6 , 690) 

690  FOHHAT('0',8  (•-')  ,  •  KG*N**2  •  ,  9  ( •  -  • )  ,/) 

WRITE  (6,700)  (WBCPMI(K)  ,K=1,3)  ,  ({WBCPAX  (K,L)  ,K=  1 , 3)  ,  L»1 , 3)  , 

6  PERFI, IERROR 

700  FORMAT  ( '  •  ,3  (F7. 3 , 1 X)  , 2X, 3 (3F9. 5,  IX)  , IX, P7. 3, 2X , 13,//) 

WRIT  E  (6 , 602) 

602  FORHATC  ',130('_')) 

C 

999  CONTIHUE 
STOP 
END 

SUBROUTINE  H ATR XM ( A ,B,C, L, M, N) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  MATRXM  PERFORMS  MATRIX  MULTIPLICATION  BETWEEN  A  MATRIX 
A  (L,  M)  AND  B(M,N)  ,  PRODUCING  A  MATRIX  C(L,N). 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DIMENSION  A(l,M),B(M,N)  ,C(L,N),R(3) 

DO  40  1=1, L 
DO  2  0  J=1,N 
R (J) =0 
DO  10  K  =  1 , H 

10  R(J)=A(I,K)  *B(K,J)  ♦R(vT) 

20  CONTINUE 
DO  30  J= 1 ,N 
30  C(I,J)*R(J) 

40  CONTINUE 
RETURN 
END 

FUNCTION  VMAG  (V) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC;  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
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C  FUNCTION  7 NAG  HEADS  A  7ECT0R  7  AND  COMPUTES  ITS  MAGNITUDE. 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

DIMENSION  7(3) 

7MAG  =  SQBT  (7  { 1)  **2*7  (2)**2*V  (3)*  *2) 

BETORN 

END 

SUBROUTINE  UNI7EC  (7,U7) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCOXCCCCCCCCCCCC 

c 

C  SUBBOOTTNB  0NI7EC  BEADS  A  7ECTOB  7  AND  COMPUTES  THE  UNIT  7 ECTOR  07 
C  ALONG  THE  SANE  DIRECTION  AS  7. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DIMENSION  7(3)  ,UV  (3) 

P  1  =  7 HAG  (7) 

IF(P1.EQ.O.)  GO  TO  20 
DO  10  K=1,3 
07  (K)=7(K)/P1 
10  CONTINUE 
GO  TO  40 
20  DO  30  K=1,j 
07 (K) =0. 

30  CONTINUE 
40  BETOBN 
END 

SUBROUTINE  7ECTP (7 1 ,72 , 73} 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  SOBBOUTINE  7ECTP  COMPOTES  THE  7ECTOR  (CROSS)  PRODUCT  (73)  OP  TNO 
C  7 ECTOR S  (71  AND  72)  . 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DIMENSION  71(3)  ,V2(3)  ,73(3) 

73(1) -(71  (2)  *72(3)) -(71  (3)  *72  (2)) 

73  (2)  =  (71  (3)  *72  (1)  )-(7l  (1)  *72(3)) 

73  (3)  =  (71(1)  *72  (2)  )  -  (7 1  (2)  *72  ( 1)  ) 

BETOBN 

END 

SUBROUTINE  TRANSP (A, AT) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  SUBROUTINE  TRANSP  TAKES  A  3  BY  3  MATRIX  (A)  AND  DETERMINES  ITS 
C  TRANSPOSE  (AT). 

C 

cCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DIMENSION  A  (3,3)  ,AT(3,3) 

DO  20  1=1,3 
DO  10  J= 1 ,3 
10  AT(I,J)=A(J,I) 

20  CONTINUE 
BETOBN 
END 

SUBROUTINE  PKDAT  (NPACK,NSEX) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  SUBROUTINE  PKDAT  READS  IN  7ARIOOS  DATA  POB  THE  4  PACKS 
C  AND  CALCULATES  CG  LOCATION  RSLATITE  TO  THE  EXTENAL  XIZ 
C  COORDINATE  SYSTEM. 

C 
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CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC'JCCCCCCCCCCCCCCCCCCCCX'CC 
COB BOB  /PK/  PKHI  (4,3)  ,PKHASS  (4)  ,  PKDIN  (4,3)  ,PXPROP(4,3)  , 

S  PK  PRO  D  (4,  3) 

COHHON  /INER/  S EG BIL ( 1 2)  , SEGB IT ( 1 2) , S EGKT D ( 1 2) , SHASS ( 12)  , 

S  PRODIN  (12,3) 

COHHON  /COFHF/  CF(4,3) 

CONBON  /COOR/  COO  (2,  12,3)  ,SHOUL  (2,3)  ,HIPPT1  (3)  ,RirPT2  (3)  ,RIFPT3<3) 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

c  READ  IN  PACK  HASSES,  HOBEKTS  AND  PRODUCTS  OP  INERTIA,  MAXIMUM 
C  D IBENS IONS  IN  THREE  DIRECTIONS,  PACK.  LOCATION  RELATIV  E  TO  SHOUL- 

C  DEH  (CORK  FACTOR) ,  AND  CG  LOCATIONS  REPRESENTED  AS  A  PERCENTAGE 
C  OF  THESE  DIMENSIONS.  THESE  PERCENTAGES  ARE  REPRESENTED 

C  IN  THE  FOLLOWING  WAY :  X-DIB-DISTANCE  IS  MEASURED  FROH 

C  FROH  THE  FRONT  OF  THB  PACK  (NEAREST  THE  BODY)  ,I-DIR- 

C  DISTANCE  IS  FROH  LEFT  PACK  EDGE,  Z-DIR:  DISTANCE  13  FROM 

C  TOP  OF  THE  PACK. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
RBAD(5,  100)  (PKHASS  (I)  ,1=1 ,4) 

100  FORHAT  (4F5.2) 

READ  (5,  110)  ((PKHI  (I,  J)  ,J=  1,3)  ,1=1,4) 

110  FORMAT  ( 3F5. 3, 3 (/, 3F5. 3) ) 

READ (5,120)  ( (P  KPROD (I , J) , J= 1 , 3)  , 1= 1 , 4) 

120  FORM  AT  ( 3 F7. 4, 3 (/, 3F7. 4) ) 

READ  (5,110)  ((PKDIH(I,  J)  ,J=1 , 3)  ,1=1,41 
READ  (5,130)  ( (P KPROP ( I  ,  J)  ,  J=  1 , 3)  ,1=1 , 4) 

130  FORMAT  ( 3F4.  2 , 3  (/ ,  3  F4.  2)  ) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  CALCULATE  CG  LOCATION  OF  SPECIFIED  PACK  RELATIVE  TO  THE 
C  EXTERNAL  XYZ  COORDINATE  ST3TEH. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
COO(NSEX, 10, 1)  =  -CP(NPACK, 1) - (PKPROP  (NPACK, 1)  *  PXDIH  (NPACK, 1) ) 
COO  ( BS6X,  10,  2)  =  (PKDi'N  (NPACK, 2)  /2.)  -  (PKDlH  (NPACK,2)  * 

6  PKPROP(NPACK,2)  ) 

COO  (NSEX,10,3)  =  GHOUL  (1,3)  *CF  (NPACK,3)  -  (PKDIS(NPACK,3)  * 

C  PK PR OP (NPACK, 3) ) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  ASSIGN  VALUES  FOR  HASS  AND  MOMENT  OF  INERTIA  FOB  SELECTED 
C  PACK  TO  ARRAYS  SHASS,  SEGHIL,  ETC.  AND  RETURN 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
SHASS  (10)  =  PKHASS  (NPACK) 

SEGHIL  (10)  =  PKHI  (NPACK,  1) 

S EGB  IT  (10)  =  PKHI  (NPACK,  2) 

SEGHID(IO)  =  PKHI  (NPACK, 3) 

DO  140  1=1,3 

PRODIN  (10,1) =PKPROD (NPACK, I) 

140  CONTINUE 
RETURN 
END 

SUBROUTINE  LDDAT  (NPACK, NSEX) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  SUBROUTINE  LDDAT  READS  IN  VALUES  FOR  THE  MAGNITUDE  OF  THE 
C  EXTRA  LOADS  AND  CALCULATES  THE  LOCATION  OF  THE  EXTRA  LOADS 

C  RELATIVE  TO  THE  EXTERNAL  XYZ  COORDINATE  SYSTEM. 
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c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

CONN  ON  /PK/  PKAI(4,3)  ,PKRASS{4)  , PKDI B (4 , 3 ) , PK PROP (4 ,3) , 

6  PKPROD  (4,  3) 

COSRON  /IVES/  SEGHIL(12) ,SEGRTT ( 12) , SEGRID( 12) ,  SHASS(12) , 

6  PRODIN  (12,3) 

CORSON  /COPRF/  CP  (4,3) 

CORSON  /COOR/  COO  (2,12,3)  ,SHOUL  (2,3)  ,RIPPT1  (3)  ,RIFPT2(3)  ,PIPPT3(3) 
REAL  LDPROP (2,3) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c  READ  IN  HAGNITODES  OP  DESIRED  EXTRA  LOADS  AND  LOCATION 

C  AS  A  PROPORTION  OF  PACK  DISENSIONS. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
READ(5,  100)  (SHASS  (I)  ,1  =  11,  12) 

100  FORR AT (2P5. 2) 

READ  (5,  110)  ( (LDPROP  (I  ,  J)  ,3=1 ,3)  ,1  =  1,2) 

110  FORR AT  (3P4. 2,/,3F4.2) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  INITIALIZE  VARIABLES 

c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  115  1=11,12 
DO  115  J  =  1 , 3 
COO (NS EX, I, J) =0. 


115  CONTINUE 

SEGRIL  (11)  =  0. 
SEGRIT(II)  =  0. 
SEGBID  (  1 1)  =  0. 
SEGRIL  (12)  =  0. 
SEGA  IT  (12)  =  0. 
SBGRID  (12)  =  0. 


DO  120  1*1,3 
PRODIN(11,I)  =0. 

PROD IN  (12, 1) =0. 

120  CONTINUE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  CALCULATE  LOCATION  OF  EXTRA  LOADS  RELATIVE  TO  THE  EXTERNAL 
C  XTZ  COORDINATE  SISTER  AND  THEN  RETURN. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
IF (SR ASS  (11) .EQ.O)  GO  TO  130 

COO(NSEX,  11,  1)  =  -CF(NPACK,1)-(LDPROP(1,1)*PKDIS(NPACK,1)) 

COO (N SEX, 11, 2)  =  (PKDIH(NPACK,2)/2.)-(PKDIR  (N PACK, 2)  * 

G  LDPROP  (1,2)) 

COO  (NSEX,  11,3)  =  SHOUL  (1,3)  tCF  ( NPACK  ,  3)  —  (PKOIH  (SPACE, 3)  * 

G  LD  PROP  (1,3)) 

130  CONTINUE 

IF(SHASS(12) .EQ.O)  GO  TO  140 

COO  (NSEX,  12,  1)  =  -CF(NPACK,1)-(LDPR0P(2,1)*PKDIH(fPACK,1>) 

COO  (NSEX, 12 ,2)  =  ( PKDIH  ( NPACK  , 2) /2.)  —  (PKDI A  (NPACK, 2)  * 

C  LDPROP  (2,2)) 

COO(NSEX,  12,  3)  =  SHOUL  (1,3)  *CF  (NPACK  ,  3)  -  (PKDI H (NPACK,  3)  * 

G  LDPROP  (2,3)) 

140  CONTINUE 
RETURN 
END 
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SUBROUTINE  AXLOC ( A , B, C , AI , A J , A K , OP) 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  SUBROUTINE  TO  LOCATE  A  UNIT  VECTOR  TBIAD  IN  A  RIGID  BODY  GIVEN 
C  THE  COORDINATES  OP  THBEE  POINTS  A,B  S  C  WHICH  ABE  NON-COLIHEAR 
C 

C  X  AXIS  HILL  BE  ALONG  AB 

C  T  AXIS  BILL  BE  THROUGH  C  PER  PEN DICUL AB  TO  AB 

C  Z  AXIS  BILL  BE  {  I  X  J  ) 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DIHENSION  A  ( 3)  ,B(3),C(3),AI(3),AJ(3)  ,  AK  (3 )  ,  OP  (3)  , TV (3) 

*  fAB(3),AC(3) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  1. DEFINE  VECTORS  AB  AND  AC 

c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL  VCOHP ( AyB, AB) 

CALL  VCOHP (A#C, AC) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  2.PIND  THE  ANGLE  BETWEEN  AB  AND  AC 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
THETA=ANGLE(AB,  AC) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  ’.CALCULATE  THE  LENGTH  OF  VECTOR  AO* 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

ACLEN=VHAG(AC) 

AOP=ACLEN*COS (THETA) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  4. ESTABLISH  THE  COORDS  OP  O* 

c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
ABI.EN=VHAG  (AB) 

PACT=  AOP/ABLEN 

CALL  INTPT(A,B,OP, FACT) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCtCCCCCCCCCCCCCCC 

c 

C  5. CALCULATE  THE  UNIT  VECTORS  ALONG  O'B  6  O'C  I 

c  I 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL  UNI VBC ( AB#  AI) 

CALL  VCOHP(OP,C,TV) 

CALL  UNI VEC (TV , A J) 


C  6. CALCULATE  UNIT  VECTOR  K  BY  VECTOB  PRODUCT 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL  VECTP(AI, AJ,AK) 

RETURN 

END 

SUBROUTINE  INTPT ( A , B, C, HAT) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 
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C  SUBROUTINE  TO  FIND  THE  COORDINATES  OP  POINT  C,  INTERMEDIATE 
C  BETWEEN  POINTS  A  AND  B  StJCH  THAT  THE  RATIO  AC:  AB  =  RAT 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

DIMENSION  A  (3)  ,B(3),C(3) 

c 

C(  1)  =  A  ( 1)  ♦  (B(1)-A{1)  )  *RAT 
C  (2)  =  A  (2)  ♦  (B{2) -A  (2)  )  *RAT 
C  (3)  =A  (3)  ♦  (B  (3 )  -  A  (3)  )  *RAT 
C 

RETORN 

END 

SUBROUTINE  VC01P(A,B,V) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  SUBROUTINE  TO  FIND  THE  COMPONENTS  OF  THE  VECTOR  JOINING  POINT  A  G  B 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DIMENSION  A  (3)  ,  B(3)  ,V  (3) 

V(1)  =  B  (  1)  -A  ( 1) 

V  (2)  =B  (2) -A  (2) 

V  (3)  =B  (3)  -A  (  3) 

RETURN 

END 

FUNCTION  ANGLE (A, B) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  PONCION  TO  FIND  THE  ANGLE  BETWEEN  TWO  VECTORS  A  AND  B. 

c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DIME NSIOR  A (3)  fB(3) ,X (3) 

V A=V H AG  (A) 

VB=VM  AG  (B) 

IF  (VA.NE.O.  .  AND.  VB.  NE.  0.  )  GO  TO  50 
C 

WRITE  (5,  10)  VA,  VB 

10  FORMAT  ( *  ***ERROR ***  IN  FUNCTION  ANGLE — DIVIDE  BY  ZERO*/ 

*  •  VMAG  (A)  =  '  ,  F 10. 3, •  VMAG(B)  =  »,F10.3) 

ANGL  B=0. 

RETURN 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  THERE  ABE  TWO  WAYS  TO  CALCULATE  THE  ANGLE  BETWEEN  TWO  VECTORS. 

C  WHEN  THE  ANGLE  IS  CLOSE  TO  0  OR  180  DEGREES,  ARSIN  IS  MORE  ACCURATE. 

C  WHEN  THE  ANGLE  IS  CLOSE  TO  90  DEGREES,  ABCOS  IS  MORE  ACCURATE. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
50  I=SC ALP  (A ,  B) 

IF (Y.GT.0.5) GO  TO  70 
ANGLE=ARCOS(Y/(VA*VB)  ) 

RETURN 

C 

70  CALI.  VECTP  ( A  ,B  ,  X) 

.  AIIGL  E=ABSIH  (VMAG  (X)  /  (VA*VB)  ) 

RETURN 

END 

FUNCTION  SCALP  (D,  E) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  FUNCTION  SCALP  COMPUTES  THE  SCALAR  (DOT)  PRODUCT  OF 
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C  TWO  VECTORS,  D  AND  F. 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

DIMENSION  D  ( 3)  ,  P(3) 

SCAL  P=D  (1)*E(1)  ♦  D(2)*E{2)*-D(3)*E(3) 

RETURN 

END 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  THE  FOLLOWING  ARE  SOURCE  LISTINGS  OF  SIX  INSL  SUBROUTINES 
C  REQUIRED  IN  THIS  PHOGRAH.  THESE  WERE  NOT  HFITTEN  BY  THE 

C  AUTHORS  BUT  ARE  CANNED  PROGFAMS  AVAILABLE  TO  PENN  ST1TB 

C  COMPUTATION  CENTER  USERS. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
IHSL  ROUTINE  NAME  -  EIGRS 


COMPUTER  -  IBM/DOUBLE 

LATEST  REVISION  -  JUNE  1,  1980 

PURPOSE  -  EIGENVALUES  AND  (OPTIONALLY)  EIGENVECTORS  OF 

A  REAL  SYMMETRIC  MATRIX 


USAGE 

ARGUMENTS  A 


N 

JOHN 


D 

Z 


IZ 


-  CALL  EIGRS  ( A,  N,  JOBN, D, Z, IZ, WK , IBR) 

-  INPUT  REAL  SYMMETRIC  MATRIX  OF  ORDER  N, 

WHOSE  EIGENVALUES  AND  EIGENVECTORS 
ABE  TO  BE  COMPUTED.  INPUT  A  IS 
DESTROYED  If  IJOB  IS  EQUAL  TO  0  OR  1. 

-  INPUT  ORDER  OF  THE  MATRIX  A. 

-  INPUT  OPTION  PARAMETER.  IP  JOBN. GE. 10 
A  IS  ASSUMED  TO  BE  IN  FULL  STORAGE  MODE 

(IN  THIS  CASE,  A  MUST  BE  DIMENSIONED  EXACTLI 
N  BY  N  IN  THE  CALLING  PROGRAM). 

IF  JOBN. LT. 10  THEN  A  IS  ASSUMED  TO  BE  IN 
SYMMETRIC  STORAGE  MODE.  DEFINE 
I JOB=HOD ( JOBN, 10) .  THEN  WHEN 

IJOB  =  0,  COMPUTE  EIGENVALUES  ONLY 
IJOB  =  1,  COMPUTE  EIGENVALUES  AND  EIGEN¬ 
VECTORS. 

IJOB  =  2,  COMPUTE  EIGENVALUES,  EIGENVECTORS 
AND  PERFORMANCE  INDEX. 

IJOB  -  3,  COMPUTE  PERFORMANCZ  INDEX  ONLY. 

IP  THE  PERFORMANCE  INDEX  IS  COMPUTED,  IT  IS 
RETURNED  IN  WK(1).  THE  ROUTINES  HAVE 
PERFORMED  (WELL,  SATISFACTORILY,  POORLY)  IF 
KK(1)  IS  (LESS  THAN  1,  EETWFEN  1  AND  100, 
GEEATER  THAN  100) . 

-  OUTPUT  VECTOR  OF  LENGTH  N, 

CONTAINING  THE  EIGENVALUE’S  OF  A. 

-  OTITPUT  N  BY  N  MATRIX  CONTAINING 

THE  EIGENVECTORS  OF  A. 

THE  EIGENVECTOR  IN  COLUMN  J  OF  Z  CORRES¬ 
PONDS  TO  THE  EIGENVALUE  D ( J) . 

IF  IJOB  -  0,  Z  IS  NOT  USED. 

-  INPUT  ROW  DIMENSION  OF  MATRIX  Z  EXACTLI  AS 

SPECIFIED  IN  THE  DIMENSION  STATEMENT  IN  THE 
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CALLING  PROGRAM. 

UK  -  WORK  AREA,  TUF  LENGTH  OF  WK  DEPENDS 

ON  THE  VALUE  OF  IJOB,  WHEN 

IJOB  =  0,  THE  LENGTH  OP  WK  TS  AT  LEAST  N. 

IJOB  =  1,  THE  LENGTH  OP  WK  IS  AT  LEAST  N. 

IJOB  =  2,  THE  LENGTH  OP  WK  IS  AT  LEAST 

N(N*1)/2*N. 

IJOB  *  3,  THE  LENGTH  OP  WK  IS  AT  LEAST  1. 

IER  -  EBBOR  PARAMETER  (OUTPUT) 

TERMINAL  ERROR 

IBB  =  128+J,  INDICATES  THAT  EQRT2S  FAILED 
TO  CONVERGE  ON  EIGENVALUE  J.  EIGENVALUES 
AND  EIGENVECTORS  1,...,J-1  HAVE  BEEN 
COMPUTED  CORRECTLY,  BUT  THE  EIGENVALUES 
ABE  UNOBDERED.  THE  PERFORMANCE  INDEX 
IS  SET  TO  1000.0 
HARNING  ERROR  (WITH  FIX) 

IN  THE  FOLLOWING,  IJOB  =  HOD (JOBN, 10) . 

IBB  =  66,  INDICATES  IJOB  IS  LESS  THAN  0  OB 
IJOB  IS  GREATER  THAN  3.  IJOB  SET  TO  1. 

IER  =  67,  INDICATES  IJOB  IS  MOT  EQUAL  TO 
ZERO,  AND  IZ  IS  LESS  THAN  THB  ORDER  OP 
MATRIX  A.  IJOB  IS  SET  TO  ZERO. 

PRECISION/ HARDWARE  -  SINGLE  AND  DOUBL E/H32 

-  SINGLE/H36, H48,H60 

REQD.  IHSL  ROUTINES  -  EHOBKS, EHOU SS , EQRT2S , UERTST , UGETIO 

NOTATION  -  INFORMATION  ON  SPECIAL  NOTATION  AND 

CONVENTIONS  IS  AVAILABLE  IN  THB  MANUAL 
INTRODUCTION  OB  THROUGH  IHSL  ROUTINE  UHELP 

COPTRIGHT  -  1980  BI  IMSL,  INC.  ALL  RIGHTS  RESERVED. 

WARRANTY  -  IHSL  WARRANTS  ONLI  THAT  IHSL  TESTING  HAS  BEEN 

APPLIED  TO  THIS  CODE.  NO  OTHER  WARRANTY, 
EXPRESSED  OR  IMPLIED,  IS  APPLICABLE. 


C 

C 

C 

C 

C 


SUBROUTINE  EIGRS 
INTEGER 

DOUBLE  PRECISION 

INTEGER 

INTEGER 

DOUBLE  PRECISION 

1 

DATA 
DATA 
♦  0/ 


IER  *  0 
JER  *  0 

IP  ( JOBN.LT.  10)  GO 
K  =  1 


(A,N, JOBN, D,Z,IZ,W K, IER) 

SPECIFICATIONS  FOR  ARGUMENTS 
V, JOBN,IZ, IER 
MU  #D(1)  , WK  (1)  ,Z  (IZ,1) 

SPECIFICATIONS  FOR  LOCAL  VARIABLES 
IJOB ,IR,JR,IJ,JI, NP1 
JER ,NA,MD,IIZ,IBEG, XL, KK  , LK, I, J, K, L 
ANORH, ASUH, PI, SUHZ,SOHR, AN,S,TEN,RDELP,ZERO, 
ONE,THOUS 

BDELP/Z341 0000000000000/ 

ZERO, ONE/O. 0D0, 1. 0D0/, TEN/10. 0D0/,THOUS/1000.0 

INITIALIZE  ERROR  PARAMETERS 
FIRST  EXECUTABLE  STATEMENT 


TO  15 

CONVERT  TO  SYMMETRIC  STORAGE  MODE 
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nno  non  no 


c 


c 


c 


c 


JI  =  N-1 
IJ  =  1 
DO  10  J*1,N 
DO  5  I=1#J 

i(K)  =  A(IJ) 

IJ  =  IJ»1 
K  =  K  ♦  1 
5  CONTINUE 

IJ  *  IJ  ♦  JI 
JI  =  JI  -  1 
10  CONTINUE 
15  IJOB  =  HOD ( JOBN, 10) 

IF  (IJOB.GE. 0. AND. IJOB. LE. 3)  GO  TO  20 

WARNING  ERROR  -  IJOB  IS  HOT  IN  THU 
BARGE 

IPR  =  66 
IJOB  =  1 
GO  TO  25 

20  IF  (IJOB. EQ. 0)  GO  TO  35 
25  IF  (IZ.GE.  N)  GO  TO  30 

WARNING  ERROR  -  IZ  IS  LESS  THi  N  N 
EIGENVECTORS  CAN  NOT  BE  COMPUTED, 
IJOB  SET  TO  ZERO 

IER  =  67 
IJOB  =  O 

30  IF  (IJOB.  EQ.  3)  GO  TO  75 
35  NA  =  (N*(N*1))/2 

IF  (IJOB. HE.  2)  GO  TO  45 
DO  40  1= 1 , N A 
WK(I)  *  A  (I) 

40  CONTINUE 

SAVE  INPUT  A  IF  IJOB  =  2 

45  ND  =  1 

IF  (I JOB. EQ. 2)  ND  =  NA+1 

REDUCE  A  TO  SIHHETBIC  TRIDI AGONAL 
FOfiH 

CALL  EHOUSS  (A,  N,  D,WK  (  .3D)  ,HK  (ND)  ) 

IIZ  -  1 

IF  (IJOB.  EQ.  0)  GO  TO  60 
IIZ  =  IZ 

SET  Z  TO  THE  IDENTITY  HATRIX 

DO  55  1=1, N 
DO  50  J=1,N 

Z(I,J)  =  ZERO 
50  CONTINUE 

Z  (1,1)  =  ONE 
55  CONTINUE 

CONFUTE  EIGENVALUES  AND  EIGENVECTORS 
60  CALL  EQRT2S  (D, WK ( ND) , N, Z, I IZ , J ER) 

IF  (IJOB. EQ. 0)  GO  TO  9000 
IF  (JER.GT.  128)  GO  TO  65 

BACK  TRANSFORM  EIGENVECTORS 

CALL  EHOBKS  (A  ,  N,  1  ,  N,  Z  ,  IZ) 

65  IF  (IJOB.LE.  1)  GO  TO  9000 

HOVE  INPUT  HATRIX  RACK  TO  A 

DO  70  1=1, NA 
A  (I)  =  WK(I) 

70  CONTINUE 

BK ( 1 )  =  THOUS 
IF  (JER.NE.  0)  GO  TO  9000 
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COMPUTE  1  -  BORB  OF  A 

75  ABORH  »  ZERO 
IBSG  *  1 
DO  85  1-1, B 
ASOfl  -  ZERO 
IL  -  IBEG 
KK  -  1 
DO  80  L*  1 ,  H 

ASOM  =  ASn.1*DA8S  (A  (IL) ) 

IF  (L.GE.I)  KK  *  l 
II  *  IL+KK 
80  COBTIBOE 

AVOBB  *  DHAX 1 ( AB03H , ASOfl) 

IBEG  =  IBEG+I 
85  COBTIBOE 

IF  f ARORB.EQ.ZERO)  ABOBB  =  OBE 

COMPOTE  PEBFORMABCE  IBDBI 

PI  *  ZERO 
DO  100  1*1, B 

IBEG  =  1  f 

S  *  ZERO 
S0RZ  *  ZEBO 
DO  95  1*1, B 
LK  *  IBEG 
KK  *  1 

SUHZ  *  SrJBZ*  DABS  (Z  (I,  I)  ) 

SOAR  *  -D  Cl) 

DO  90  K=1 ,  B 

SUM B  *  SHSRt  A  (LK)  *Z  (K,  I) 

IF  (K.GE.L)  KK  *  K 

IK  =  IK*  KK  j 

90  COBTIBOE 

S  =  S*CABS  (SOAR) 

IBEG  =  IBEG  *L 
95  COBTIBOE 

IF  (SOBZ.EQ.ZERO)  GO  TO  100 
PI  *  DBAK1 (PI,S/SOnZ) 

100  COBTIBOE 
AB  *  B 

PI  *  PI/(AHORH*TSB*AM*RDELP) 

WK(1)  =  PI 

IF  (JOBB.LT. 10)  GO  TO  9000 

COBVEBT  BACK  TO  FOLI  STORAGE  BODE 

bpi  *  n*l 

IJ  =  (B-l)^BPl  ♦  2 
K  *  (B*(BP1)  )/2 
DO  110  JR* 1 , B 
J  =  BP1-JR 
DO  105  IR*1,J 
IJ  *  I J-1 
A  (IJ)  *  A  (K) 

K  *  K-1 

105  COBTIBOE 

IJ  *  IJ-JR 
110  COBTIBUE 
JI  =  0 
K  *  M-1 
DO  120  1*1,  B 
IJ  *  I-B 
DO  115  J*1,I 
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IJ  =  IJ  +  N 
JI  =  JIM 
A  (IJ)  =  A  (JI) 

115  CONTINUE 

JI  =  JI  ♦  K 
K  =  K-1 
120  CONTINOE 
9000  CONTINUE 

IP  (IER.NE.0)  CALL  UEETST  (IER,6HSIGRS  ) 
IP  (JER.EQ.O)  GO  TO  9005 
IER  =  JER 

CALL  OEHTST  (I  ER ,  6  HEIGRS  ) 

9005  RETURN 
END 

IMSL  ROUTINE  NAME  -  EHOBKS 


COMPUTER 
LATEST  REVISION 
PURPOSE 

USAGE 

ARGUMENTS  A 

N 

Ml 


M2 

Z 


IZ 


PRECISION/HARDWARE 


-  IBM/DOUBLE 

-  JANUARY  1,  1978 

-  BACK  TRANSFORMATION  TO  FORM  THE  EIGENVECTORS 

OF  THE  ORIGINAL  SYMMETRIC  MATRIX  FROM  THE 
EIGENVECTORS  OF  THE  TRIDTAGONAL  MATRIX 

-  CALL  EHOBKS  ( A, N , H 1 , M2, Z, IZ) 

-  THE  AHRAY  CONTAINS  THE  BETA  TLS  OF  THE  HOUSE¬ 

HOLDER  REDUCTION  OF  THE  ORIGINAL  MATRIX  A 
AS  GENERATED  BY  IMSL  ROUTINE  EIIOUSS.  (INPUT) 

-  ORDER  OF  THE  REAL  STHMETRIC  MATRIX.  (INPUT) 

-  Ml  AND  M2  ARE  TWO  IUPUT  SCALARS  SUCH  THAT 

EIGENVECTORS  Ml  TO  M2  OF  THE  TRIDIAGONAL 
MATRIX  A  HAVE  BEEN  FOUND  AND  NORMALIZED 
ACCORDING  TO  THE  EUCLIDEAN  HORN. 

-  SEE  ABOVE  -  Ml 

-  A  TWO  DIMENSIONAL  ARRAY  OF  SIZE  N  X  (K2-M1  +  1) 

WHICH  CONTAINS  EIGENVECTORS  HI  TO  M2  OF 
TRIDI AGON  AL  MATRIX  T,  NORMALIZED  ACCORDING 
TO  EUCLIDEAN  HORN.  INPUT  Z  CAN  BE  PRODUCED 
BY  IMSL  ROUTINE  EQR  T2Sr  THE  RESULTANT 
MATRIX  OVERWRITES  THE  INPUT  Z.  (INPUT/OUTPUT) 

-  ROW  DIMENSION  OF  MATRIX  Z  EXACTLY  AS 

SPECIFIED  IN  THE  DIMENSION  STATEMENT  IN  THE 
CALLING  PROGRAM.  (INPUT) 

-  SINGLE  AND  DOUBLE/H32 

-  SIHGLE/H36, H48,H60 


R  EQD.  IMSL  ROUTINES  -  NONE  REQUIRED 

NOTATION  -  INFORMATION  ON  SPECIAL  NOTATION  AND 

CONVENTIONS  IS  AVAILABLE  IN  THE  MANUAL 
INTRODUCTION  OR  THROUGH  IMSL  ROUTINE  UHELP 

COPYRIGHT  -  1978  EY  IMSL,  INC.  ALL  RIGHTS  RESERVED. 

WARRANTY  -  IMSL  WARRANTS  ONLY  THAT  IHSL  TESTING  HAS  BEEN 

APPLIED  TO  THIS  CODE.  NO  OTHER  WARRANTY, 
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EXPRESSED  OB  XHPLIBD,  IS  APPLICABLE 


C 
C 
c- 
c 

SUBROUTINE  EHODKS  (A, N  ,  H 1 , 82 , Z , IZ) 
C 

DIME NSIOM  A  (1)  ,Z  (IZ,  1) 

DOUBLE  PRECISION  A,Z,H,S 
C 

IP  (N  .BQ.  1)  GO  TO  30 

DO  25  1*2, N 
L  *  1-1 
IA  *  (I*L)/2 
H  *  A(IA*I) 

IP  (H.EQ.0.D0)  GO  TO  25 
C 
C 
C 
C 

DO  20  J  *  HI, M2 
S  ■  O.ODO 
DO  10  K  *  1 , L 

S  *  S*A(IA>K)  *Z  (K,J) 

10  CONTINUE 

S  *  S/H 
DO  15  K-1,L 

Z  (K ,  J)  *  Z  (K,  J)  -S*A(IA*K) 
15  CONTINUE 

20  COMTIIUF 
25  CONTINUE 
30  RETURN 
END 

IHSL  ROUTINE  NAME  -  EHOUSS 


FIRST  EXECUTABLE  STATEMENT 


DERIVES  EIGENVECTORS  81  TO  82  OF 
THE  ORIGINAL  MATRIX  PROM  EIGENVECTORS 
HI  TO  82  OP  THE  SIHHETHIC 
TBIDIAGORAL  MATRIX 


COMPUTER  -  IBM/DOUBLE 

LATEST  REVISION  -  JANUARY  1,  1978 

PURPOSE  -  REDUCTION  OP  A  SYMMETRIC  MATRIX  TO  SYMMETRIC 

TRIDIAGONAL  FORM  USING  A  HOUSEHOLDER 
DEDUCTION 


OSAGE 


-  CALL  EHOUSS  (A,N,D,E,E2) 


ARGUMENTS  A 


N 

D 


P.2 


-  THE  GIVEN  N  X  N,  PEAL  SYMMETRIC  MATRIX  A, 

WHERE  A  IS  STORED  IN  SYMMETRIC  STORAGE  MODE. 
THE  INPUT  A  IS  REPLACED  BY  THE  DETAILS  OP 
THE  HOUSEHOLDER  REDUCTION  OP  A. 

-  INPUT  ORDER  OP  A  AND  THE  LENGTH  OP  D,  E,  AND 

E2. 

-  THE  OUTPUT  ARRAY  OF  LENGTH  N,  GIVING  THE 

DIAGONAL  ELEMENTS  OP  THE  TRIDIAGONAL  MATRIX. 

-  THE  OUTPUT  ARRAY  OF  LENGTH  H,  GIVING  THE  SUB- 

DIAGONAL  IN  THE  LAST  (N-1)  EL  EBERTS,  B(1)  IS 
SET  TO  ZERO. 

-  OUTPUT  ARRAI  OF  LENGTH  N.  E2  (I)  «  E|I)**2. 


P RECISION/HABDVABE  -  SINGLE  AND  DOUBLE/H32 
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-  SINGLE/H36, H48,H  60 
REQD.  IHSL  ROUTINES  -  NONE  REQU  IR ED 

NOTATION  -  INFORMATION  ON  SPECIAL  NOTATION  AND 

CONVENTIONS  IS  AVAILABLE  IN  THE  MANUAL 
INTRODUCTION  OP  THROUGH  IHSL  ROUTINE  UHELP 

COPIRIGHT  -  1978  BT  IHSL,  INC,  ALL  RIGHTS  RESBRVED. 

WARRANT!  -  IHSL  NARRANTS  ONLY  THAT  IHSL  TESTIHG  HAS  BEEN 

APPLIED  TO  THIS  CODE.  NO  OTHBB  WARRANT!, 
EXPRESSED  OR  IMPLIED,  IS  APPLICABLE. 


SUBROUTINE  BHOUSS  (A, N, D, E, E2) 

C 

DIMENSION  A(1)  ,  )(N)  ,E(N)  ,  E2  (N) 

DOUBLE  PRECISION  A, D, E, E2, Z ERO, H, SCALE, ONE, SCALE1 ,F,G, HH 
DATA  ZERO/O. 0D0/, ONE/1. 0D0/ 

C  FIRST  EXECUTABLE  STATEHENT 

NP1  =  N ♦  1 
NN  =  (N*HPl)/2-1 
NBEG  *  NH+1-N 
DO  70  II  =  1 ,N 
I  a  NP1-TI 

L  *  1-1 

H  *  ZERO 

SCALE  =  ZERO 

IF  (L  .LT.  1)  GO  TO  10 

C  SCALE  ROW  (ALGOL  TOL  THEN  NOT  NEEDED) 

NK  =  NN 
DO  5  K  =  1 ,  L 

SCALE  =  SCALE+DABS  (A  (NK) ) 

NN  =  NK- 1 
5  CONTINUE 

IF  (SCALE  .  NK.  ZERO)  GO  TO  15 
10  E  (I)  =*  ZERO 

E2  (I)  =  ZERO 
GO  TO  65 
15  NK  =  NN 

SCALE1  =  ONE/SCALE 
DO  20  K  ®  1,1 

A  (NK)  =  A  (NK)  *SCALE1 
H  =  H+A  (NK)  *  A  (NK) 

NK  =  NK-  I 
20  CONTINUE 

E2  (I)  =  SCALE*SCALE*H 
F  =  A (NN) 

G  =  -DSIGN  (DSQRT  (H)  ,  P) 

E  (I)  =  SCALB*G 

H  =  H-F*G 
A (NN)  =  F-G 

IF  (L  .EQ.  1)  GO  TO  55 
F  *  ZERO 
JK1  *  1 
DO  40  J  =  1 ,  L 
G  =  ZERO 
IK  a  NBEG+1 
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JK  -  J*1  fOBn  RLEBBNT  op  A*U 

DO  25  K  3  1  f J 

G  =  G*A(JK)*A(IK) 

JK  *  JK+1 
IK  *  IK*1 
CONTINUE 
Jpl  *  J*1 

IF  (L  .LT.  JPD  GO  T0  35 
JK  -  JK*J-1 
DO  30  K  *  JP1#I 

G  *  G*A  (JK)  *A  (IK) 

JK  *  JK+K 

IK  *  IK*1 

CONTINUE  poBR  2LEHEMT  OP  P 

E  ( J)  3  G/H 
F  «  F*B(J)  *A  (NBRG+J) 

JK1  «  JK1+J 
CONTINUE 

HR  3  r/(H*H)  ?0FH  bEDTICSO  A 


fofh  reduced  a 


JK  *  1 
DO  50  J  * 

F  s  A(NBEG*J) 

G  =  E(J)  -HH*F 

E(J)  3  G 

°°  mj£)=3  a(JK>-p*e(K)-g*a(nbeg*k) 
JK  =  JK+1 
45  CONTINUE 

50  CONTINUE 

55  00  SCiLS*S  (BBEG+K) 

go  CONTINUE 

65  !»«C»LS.SC11S 

NBBG  *  NBEG-I* 1 
NN  3  NN-T 
70  CONTINUE 
RETURN 

IHSLDROUTINE  NAHE  -  EQRT2S 


COHPUTBF 
LATEST  REVISION 
PURPOSE 


USAGE 

abguhents 


-  IBH/DOUBLE 

-  JAMUART  1,  1978 

QL  HBTHOD. 

-  CALL  BQRT2S  (D,E,N,2, IZ,IER) 

-  ON  INPUT,  THE  VECTOR  ° 

THE  DIAGONAL  ELEBEMTS  OP  THE  SIBHETHIC 

SPSSSStVSKillS  IBS  EIGEStUOES  OP 


o  n  o 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 

c 


T  IN  ASCENDING  ORDER. 

E  -  ON  INPOT,  THE  VECTOR  E  OF  LENGTH  N  CONTAINS 

THE  SUB-DIAGONAL  ELEMENTS  OF  T  IN  POSITION 
2,...,N.  ON  OOTPUT,  E  IS  DESTROTED. 

N  -  ORDER  OF  TRIDIAGONAL  MATRIX  T. (INPUT) 

Z  -  ON  INPOT,  Z  CONTAINS  THE  IDENTITY  MATRIX  OF 

ORDER  N. 

ON  OOTPOT,  Z  CONTAINS  THE  EIGENVECTORS 
OF  T.  THE  EIGENVECTOR  IN  COLUMN  J  OF  Z 
CORRESPONDS  TO  THE  EIGENVALUE  D  (.J)  . 

IZ  -  INPUT  BOH  DIMENSION  OF  MATRIX  Z  EXACTLY  AS 

SPECIFIED  IN  THE  DIMENSION  STATEMENT  IN  THE 
CALLING  PROGRAM.  IF  IZ  IS  LESS  TRAN  N,  THE 
EIGENVECTORS  ARE  NOT  COMPUTED.  IN  THIS  CASE 
Z  IS  NOT  USED. 

I  EH  -  ERROR  PARAMETER 
TERMINAL  ERROR 

IEH  =  128*J,  INDICATES  THAT  EQBT2S  FAILED 
TO  CONVERGE  ON  EIGENVALUE  J.  EIGENVALUES 
AND  EIGENVECTORS  HAVE  BEEN 

COMPUTED  CORRECTLY,  BUT  THE  EIGENVALUES 
ABE  UNORDERED. 


PRECIS  ION/HARD WARE  -  SINGLE  AND  DOUBLE/H32 

-  SINGLE/H36, H48,H60 


REQD.  IHSL  ROUTINES  -  UBitTST, UGET 10 


NOTATION  -  INFOBHSTION  ON  SPECIAL  NOTATION  AND 

CONVENTIONS  IS  AVAILABLE  IN  THE  MANUAL 
INTRODUCTION  OR  THROUGH  IHSL  ROUTINE  UHELP 


REMARKS  IHSL  POUTINE  EQRT2S  IS  DESIGNED  TO  ACCEPT  OUTPUT 

VECTORS  D  AND  E  FROM  IMSL  ROUTINE  EHOUSS  AS  INPUT 
D  AND  E  OF  EQRT2S.  GIVEN  A  SYMMETRIC  TRIDIAGONAL 
MATRIX,  T,  VECTOR  D  CONTAINS  THE  DIAGONAL  ELEMENTS 
OF  T  AND  VECTOR  E  CONTAINS  THE  SUBDIAGONAL  ELEMENTS 
OF  T.  SEE  THB  EHOUSS  DOCUMENT. 


COPYRIGHT 


1978  BY  IMSL,  INC.  ALL  FIGHTS  RESERVED. 


WARRANTY  -  IHSL  WARRANTS  ONLY  THAT  IHSL  TESTING  HAS  BEEN 

APPLIED  TO  THIS  CODE.  NO  OTHER  WARRANTY, 
EXPRESSED  OR  IMPLIED,  IS  APPLICABLE. 


SUBROUTINE  EQRT2S  (D, E , N, Z , IZ , I EF) 


DIMENSION 
DOUBLE  PRECISION 
DATA 
DATA 


D  (  1 )  ,E(1)  ,Z(IZ,  1) 

D,E,Z,B,C,F,G,H,P,R,S, RDELP, ONE, ZERO 
BDELP/Z34 100 000000 00000/ 

ZERO, ONE/O. ODO, 1.0D0/ 

HOVE  THE  LAST  N- 1  ELEMENTS 

OF  E  INTO  THE  FIRST  N-1  LOCATIONS 

FIRST  EXECUTABLE  STATEMENT 


IER  =  0 

IF  (N  .EQ.  1)  GO  TO  9005 
DO  5  I=2,N 

E(I-1)  *  E(I) 
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c 


c 


c 


5  CONTINUE 
E (M)  =  ZEBO 
B  *  Z  EB  0 
P  =  ZEBO 
DO  60  1=1, N 


J  =  0 

H  =  PDELP*(DIBS  (D(L)  )  ♦DiBS  (B  (L)  )) 
IT  tB.LT.R)  B  *  H 

LOOK  POB 


SHALL  SOB-DIAGONAL  BLBBtVT 


DO  10  E=  L,  N 

K=« 

IP  (DABS  (E  (KJ)  «I*E.  B)  GO  TO  15 
10  CONTINUE 

15  N  =  K 

IP  (B.EQ.L)  GO  TO  55 

20  IF  (J  .£Q.  30)  GO  TO  B5 

J  *  J  +  1 
Li  =  L*  1 
G  *  D(L) 

P  =  (0(L1)-G)/(E(L)+E(L)) 

F  =  DSQRT  (P*P*ONE) 

D  (L)  =  E(L)/(P+DSIGN(R,P)) 

FI  =  G-D(L) 

DO  25  I  =  L1,H 
D  (I)  *  D  (I)  -H 
25  CONTINUE 

f  * 

QL  TRANSPOSE  AT ION 


P  =  D  ((1) 

C  =  ONE 
S  =  ZERO 
(INI  =  (1-1 
NH1PL  =  NK1+L 
IP  (L.GT.  (INI )  GO  TO  50 
DO  45  II  =  L,  (1(1 1 
I  *  NH1PL-II 
G  =  C*B  (I) 

H  =  C*P 

IP  (DABS  (P)  ,LT«  DABS  (E  (I)  )  )  GO  TO  30 
C  =  B(I)/P 
B  =  DSQRT  (C*C*ONE) 

B  (1*1)  *  S*P*B 
S  =  C/R 
C  =  ONE/F 
GO  TO  35 

30  C  =  P/EfI) 

R  =  DSQRT  (u*C*OHE) 

E(I*1)  =  S*E  (I)  *F 
S  *  ONE/H 
C  *  C*S 

35  P  =  C*D(I)-S*G 

D  (1+1)  *  H*S*  (C*G*S*D  (I)  ) 

IF  (IZ  -IT.  N)  GO  TO  45 

PORN  VECTOR 

DO  40  K=1,N 

H  =  Z  (K,  1*1) 

2  (K,  I*  1)  =  S*Z{K,I)  *C*B 
Z(K,I)  =  C*Z  (K,I)  ”S*H 
40  CONTINUE 

45  CONTINUE 


79 


nonnnnnnnononononnnnnnononno 


50 


E(L)  *  S*P 
D(L)  *  C*P 

IP  (DABS  (E(L)  i  .GT.B)  GO  TO  20 
55  0  (L)  =  D(L)  ♦  F 

60  CONTI  NOE 

C  ORDER  ExGENV ALOES  AND  EIGENVECTORS 

DO  80  1=1, N 

K  =  I 
P  =  D (I) 

IP1  =  1*1 

IF  (IF1.GT.N)  GO  TO  70 
DO  65  J=IP1 , N 

IF  (D  (  J)  .  GE.  P)  GO  TO  65 
K  =  J 
P  =  D  (J) 

65  CONTINUE 

70  IP  (K.EQ.I)  vr'  TO  80 

D(IC)  *  D(I) 

D (I)  *  ? 

IF  (IZ  .LT.  N)  GO  TO  80 
DO  75  J  =  1,  N 
P  =  Z(J,I) 

Z(J,I)  =  Z(J,K) 

7,  ( J,  K)  =  P 

75  CONTINUE 

80  CONTINUE 
GO  TO  9005 
85  IER  =  1 28*L 
9000  CONTINUE 

CALL  UEBTST (IER ,6 H EQHT2S) 

9005  HETUBN 
END 

I BSL  ROUTINE  NINE  -  UEFTST 


COMPUTER 
LATEST  REVISION 
PURPOSE 
OSAGE 

A  RGUME  NTS  IER 


BANE 

PRECISION/HABDVARE 
REQD.  INSL  BOOTIWES 
NOTATION 


-  IBH/SIHGL*’ 

-  JANUARY  1,  1978 

-  PRINT  A  MESSAGE  REFLECTING  AN  ERROR  CONDITION 

-  CALL  OEBTJT  (IER  .NAME) 

-  ERROR  PARAMETER.  (INPUT) 

IER  =  I*J  WHERE 

I  =  128  IMPLIES  TERMINAL  ERRCP, 

I  =  64  IMPLIES  WARNING  WITH  FIX,  AND 

I  =  32  IMPLIES  WARRING. 

J  *  BRBOR  CODE  RELEVANT  TO  CALLING 
ROUTINE. 

-  A  SIX  CHARACTER  LITERAL  STRING  GIVING  THE 

IABE  OF  THE  CALLING  ROUTINE^  (INPUT) 

-  SI NGLE/ALL 

-  UGETIO 

-  INFORMATION  ON  SPECIAL  NOTATION  AND 

CONVENTIONS  IS  AVAILABLE  IN  THE  MANUAL 

80 


c 

c 

C  REMARKS 
C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C  COPYRIGHT 
C 

C  WARRANTY 

C 

C 

C 

C - 


I NTRODUCTION  OR  THROUGH  IHSL  IOOTIHE  UHELP 

THE  ERROR  MESSAGE  PRODUCED  BY  UFRTST  IS  WRITTEN 
ONTO  THE  STANDARD  OUTPUT  UHIT.  THE  OUTPUT  OMIT 
NUMBER  CAN  BE  DETERMINED  BY  CALLING  UGETIO  AS 
FOLLOWS..  CALL  UGETIO (1 , HIM, MOUT) . 

THE  OUTPUT  UHIT  HUMBER  CAM  BE  CHAMGED  BY  CALLING 
UGETIO  AS  FOLLOWS. . 

MIH  =  0 

HOUT  =  NEW  OUTPUT  UHIT  HUMBER 
CALL  UGETIO  (3, HIM, HOUT) 

SEE  THE  UGETIO  DOCUMENT  FOR  MORE  DETAILS. 

-  1978  BY  IHSL,  INC.  ALL  RIGHTS  RESERVED. 

>  I MSL  WARRANTS  ONLY  THAT  IHSL  TESTING  HAS  BBS! 
APPLIED  TO  THIS  CODE.  NO  OTHER  WARRANTY, 
EXPRESSED  OR  IMPLIED,  IS  APPLICABLE. 


C 

C 

C 

C 

C 

c 

C 


SUBROUTINE  HERTST  (IFR , 1AME) 

SPECIFICATIONS  FOR  ARGUMENTS 

INTEGER  IER 

IMTEGER*2  NAME  (3) 

SPECIFICATIONS  FOB  LOCAL  VARIABLES 
INTEGERS  HAHSET(3)  ,NAMEQ(3) 

DATA  NAHSET/2HUE, 2BRS, 2HET/ 

DATA  MAHEQ/2B  ,2H  ,2H  / 

FIRST  EXECUTABLE  STATEMENT 
DATA  LEVEL/4/, IEQDF/0/, IEQ/1H*/ 

IF  (IER. GT. 999)  GO  TO  25 
IF  (IEH.LT.-32)  GO  TO  55 
IF  (IER. LE.  128)  GO  TO  5 
IF  (LEVEL-LT.  1)  GO  TO  30 

PRINT  TERMINAL  MESSAGE 

CALL  UGETIO (1, HI N,IOU NIT) 

IF  (IEQDF.EQ.  1)  WHITE (IOUNIT, 35)  IER, NAHEQ, IEQ, NAME 
IF  (IEQDF.EQ.O)  HR ITE (IOUN7T ,35)  IER, NAME 
GO  TO  30 

5  IF  I  IER.  LE.  64)  GO  TO  10 
IF  ( LEVEL. LT. 2)  GO  TO  30 

PBINT  WARNING  WITH  FIX  MESSAGE 

CALL  UGETIO(1,NIN,  IOUNIT) 

IF  (IEQDF.BQ.1)  WRITE (IOUNIT, 40)  IER , NAMEQ, IEQ, NAME 
IF  (IEQDF.EQ.O)  WRITE  (IOUNIT, 40)  IER, NAME 
GO  TO  30 

10  IF  (IER.LE.  32)  GO  TO  15 

PRINT  WARNING  MESSAGE 

IF  (LEVEL. LT. 3)  GO  TO  30 
CALL  UGETIO (1,NIN, IOUNIT) 

IF  (IEQDF.EQ.1)  WRITE (IOUNIT, 45)  IER, MAMEQ,IEQ, NAME 
IF  (IEQDF.EQ.O)  WRITE ( IOUNIT , 45)  IER, NAME 
GO  TO  30 
15  CONTINUE 


CHECK  FOR  UBBSET  CALL 

DO  20  1*1,3 

IF  (NAME  (I)  .  NE.  NAM  SET  (I) )  GO  TO  25 
20  CONTINUE 

LEVOLD  *  LEVEL 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LEVEL  =  I SB 
IER  =  LEVOLD 

IF  (LEVBL.LT. 0)  LEVEL  =  4 
IF  (LEVEL. GT. 4)  LEVEL  =  4 
GO  TO  30 
25  CONTI  HOB 

IF  (LBVEL.LT. 4)  GO  TO  30 

PRINT  NON-DEFINED  MESSAGE 

CALL  UGETIO (1 ,NIN, IOONIT) 

IF  (IEQDF. EQ.1)  WHITE  (IOONIT, 50)  I RP , NA MEQ, IEQ, N A HE 
IF  (IEQDP.EQ.O)  WR ITE ( IOONIT , 50)  IER, NAME 
30  IEQDF  =  0 
PETORN 

35  FORMAT ( 19H  ***  TERMINAL  ERROR, 10X, 7H (IER  =  ,13, 

1  20H)  FROM  IHSL  BOOTINE  ,3i2,Al»3A2) 

40  FORMAT (36b  ***  WARNING  WITH  FIX  ERROR  (IER  =  ,13, 

1  20H)  FROM  IHSL  ROOTINE  ,3A2,A1,3A2) 

45  FORM  AT  ( 18H  WARNING  ERROR , 1 1 X,7H  (IER  =  ,13, 

1  20H)  FROM  IHSL  ROOTINE  ,3A2,A1,3A2) 

50  FORMAT  (20H  ***  OHDIFINED  ERROR, 9X,7H (IER  =  ,15. 

1  2 OH)  FROM  IHSL  ROUTINE  ,3A2,A1,3A2) 

SAVE  P  FOB  P  =  R  CASE 
P  IS  THE  PAGE  NAME 
R  IS  THE  BOnTINB  NAME 

55  IEQDF  =  1 
DO  60  1=1,3 
60  NAHEQ(I)  =  MANE  (I) 

65  RETOBN 
BID 

IHSL  ROUTINE  NAME  -  UGETIO 


COMPOTEF  -  IBH/SINGLE 

LATEST  REVISION  -  JANOART  1,  1978 

PORPOSE  -  TO  RETRIEVE  COBRENT  V ALOES  AND  TO  SET  NEW 

VALOES  POP  INPUT  AND  OUTPUT  ONIT 
IDENTIFIERS. 

OSAGE  -  CALL  OGET 10 (IOPT , NI N, NOUT) 

ARGUMENTS  IOPT  -  OPTION  PARAMETER.  (INPOT) 

IP  IQPT= 1 ,  THE  CURRENT  INPOT  AND  OOTPOT 
UNIT  IDENTIFIER  VALOES  ARE  RETURNED  IN  NIN 
AND  NOUT,  RESPECTIVELT. 

IF  IOPT=2  (3)  THE  INTERNAL  VALOE  OF 
NIN  (NOOT)  IS  RESET  FOR  SUBSEQUENT  USE. 

NIN  -  INPUT  UNIT  IDENTIFIER. 

OUTPUT  IF  IOPT=  1 ,  INPUT  IF  IOPT=2. 

NOUT  -  OUTPUT  UNIT  IDENTIFIER. 

OUTPUT  IF  IOPT= 1 ,  INPUT  IF  IOPT=3. 

PRECISION/ HARDWARE  -  SINGLE/ALL 
REQD.  IHSL  ROUTINES  -  NONE  REQUIRED 

NOTATION  -  INFORMATION  ON  SPECIAL  NOTATION  AND 

CONVENTIONS  IS  AVAILABLE  IN  THE  MANUAL 
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l«RODyCTlOH  OR  THROUGH  IHSL  BOOTIES  OHRT.P 
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HIED, EOUTD 

KIED/5A  «°  J«/6/Elg caT  k0LB  STITEHENT 


IF  (IOPT.  EQ.3) 
IF  (IOPT-  EQ*  2) 
IF  (XOPT.HE.  1) 
EIE  *  EIED 
ROOT  *  EOOTD 
GO  TO  9005 
5  HIED  *  «I* 

GO  TO  9005 
10  EOOTD  *  EO0T 
9005  BETOM 
t«D 


GO  TO  10 
GO  TO  5 
GO  TO  9005 


